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Abstract 



Cosmologists today are confronted with the perplexing reahty that the uni- 
verse is currently accelerating in its expansion. Nevertheless, the nature of the 
fuel that drives today's cosmic acceleration is an open and tantalizing mystery. 
There exists the intriguing possibility that the acceleration is not the man- 
ifestation of yet another mysterious ingredient in the cosmic gas tank (dark 
energy), but rather our first real lack of understanding of gravity itself, and 
even possibly a signal that there might exist dimensions beyond that which 
we can currently observe. The braneworld model of Dvali, Gabadadze and 
Porrati (DGP) is a theory where gravity is altered at immense distances by 
the excruciatingly slow leakage of gravity off our three-dimensional Universe 
and, as a modified-gravity theory, has pioneered this line of investigation. I 
review the underlying structure of DGP gravity and those phenomenological 
developments relevant to cosmologists interested in a pedagogical treatment 
of this intriguing model. 



*E-mail: lue@lonestar.utsa.edu 



I. THE CONTEMPORARY UNIVERSE 



Cosmology in this decade is said to be thriving in a golden age. With the cornucopia of 
observational data from both satellite and ground-based surveys, an increasingly coherent 
phenomenological and theoretical picture is now emerging. And while our understanding 
of cosmic expansion, primordial nucleosynthesis, the microwave background and other phe- 
nomena allows particle physics and cosmology to use the very vastness of our Universe to 
probe the most incomprehensibly high energies, this golden age of cosmology is offering the 
first new data regarding the physics on immense scales in of themselves. In other words, 
while modern cosmology is a ratification of the notion of the deep connection between the 
very small and the very large, it offers also the opportunity to challenge fundamental physics 
itself at the lowest of energies, an unexplored infrared domain. 

A central example highlighting this theme is that physicists are currently faced with the 
perplexing reality that the Universe is accelerating in its expansion [1,2]. That startling 
reality is only driven home with the observation of the onset of this acceleration [3]. The 
acceleration represents, in essence, a new imbalance in the governing gravitational equations: 
a universe filled only with ordinary matter and dark matter (ingredients for which we have 
independent corroboration) should decelerate in its expansion. What drives the acceleration 
thus remains an open and tantalizing question. 

Instructively, physics historically has addressed such imbalances in the governing gravita- 
tional equation in either one of two ways: either by identifying sources that were previously 
unaccounted for (e.g., Neptune and dark matter) or by altering the governing equations 
(e.g., general relativity). Standard cosmology has favored the first route to addressing the 
imbalance: a missing energy-momentum component. Indeed, a "conventional" explanation 
exists for the cause of that acceleration — in general relativity, vacuum energy provides 
the repulsive gravity necessary to drive accelerated cosmological expansion. Variations on 
this vacuum-energy theme, such as quintessence, promote the energy density to the poten- 
tial energy density of a dynamical field. Such additions to the roster of cosmic sources of 
energy-momentum are collectively referred to as dark energy. If it exists, this mysterious 
dark energy would constitute the majority of the energy density of the universe today. 

However, one may also entertain the alternative viewpoint. Can cosmic acceleration be 
the first new signal of a lack of understanding of gravitational interactions? I.e., is the 
cosmic acceleration the result, not of the contents of the the cosmic gas tank, as it were, but 
a consequence of the engine itself. This is the question that intrigues and excites us, and 
more importantly, how we can definitively answer that question. How can one definitively 
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differentiate this modification of the theory of gravity from dark energy? Cosmology can 
offer a fresh opportunity to uncover new fundamental physics at the most immense of scales.^ 

Understanding cosmic acceleration and whether it indicates new fundamental physics 
serves as a first concrete step in the program of exploiting cosmology as a tool for un- 
derstanding new infrared physics, i.e., physics on these immense scales. In 2000, Dvali, 
Gabadadze and Porrati (DGP) set forth a braneworld model of gravity by which our ob- 
served four- dimensional Universe resides in a larger, five- dimensional space. However, unlike 
popular braneworld theories at the time, the extra dimension featured in this theory was 
astrophysically large and fiat, rather than large compared to particle-physics scales but oth- 
erwise pathologically small compared to those we observe. In DGP braneworlds, gravity 
is modified at large (rather than at short) distances through the excruciatingly slow evap- 
oration of gravitational degrees of freedom off of the brane Universe. It was soon shown 
by Deffayet that just such a model exhibits cosmological solutions that approach empty 
universes that nevertheless accelerate themselves at late times. 

Having pioneered the paradigm of self-acceleration, DGP braneworld gravity remains a 
leading candidate for understanding gravity modified at ultralarge distances; nevertheless, 
much work remains to be done to understand its far-reaching consequences. This article 
is intended to be a coherent and instructive review of the material for those interested in 
carrying on the intriguing phenomenological work, rather than an exhaustive account of the 
rather dissonant, confusing and sometimes mistaken literature. In particular, we focus on the 
simple cases and scenarios that best illuminate the pertinent properties of DGP gravity as 
well as those of highest observational relevance, rather than enumerate the many intriguing 
variations which may be played out in this theory. 



There is a bit of a semantic point about what one means by dark energy versus modified gravity, 
i.e., altering the energy-momentum content of a theory versus altering the field equations them- 
selves. For our qualitative discussion here, what I mean by dark energy is some (possibly new) field 
or particle that is minimally coupled to the metric, meaning that its constituents follow geodesies of 
the metric. An alternative statement of this condition is that the new field is covariantly conserved 
in the background of the metric. I am presuming that the metric alone mediates the gravitational 
interaction. Thus, a modified-gravity theory would still be a metric theory, minimally coupled to 
whatever energy-momentum exists in that paradigm, but whose governing equations are not the 
Einstein equations. 

We also wish to emphasize the point that we are not addressing the cosmological constant problem 
here, i.e., why the vacuum energy is zero, or at least much smaller than the fundamental Planck 
scale, Mp. When we refer to either dark energy or a modified-gravity explanation of cosmic 
acceleration, we do so with the understanding that a vanishing vacuum energy is explained by 
some other means, and that dark energy refers to whatever residual vacuum energy or potential 
energy of a field may be driving today's acceleration, and that modified-gravity assumes a strictly 
zero vacuum energy. 
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We begin by setting out the governing equations and the environment in which this 
model exists, while giving a broad picture of how its gross features manifest themselves. We 
then provide a detailed view of how cosmology arises in this model, including the emergence 
of the celebrated self-accelerating phase. At the same time, a geometric picture of how such 
cosmologies evolve is presented, from the perspective of both observers in our Universe as 
well as hypothetical observers existing in the larger bulk space. We touch on observational 
constraints for this specific cosmology. We then address an important problem/subtlety 
regarding the recovery of four-dimensional Einstein gravity. It is this peculiar story that 
leads to powerful and accessible observable consequences for this theory, and is the key to 
differentiating a modified-gravity scenario such as DGP gravity from dark-energy scenarios. 
We then illuminate the interplay of cosmology with the modification of the gravitational 
potentials and spend the next several sections discussing DGP gravity's astronomical and 
cosmological consequences. Finally, we finish with the prospects of future work and some 
potential problems with DGP gravity. We will see that DGP gravity provides a rich and 
unique environment for potentially shedding light on new cosmology and physics. 

II. GRAVITATIONAL LEAKAGE INTO EXTRA DIMENSIONS 

We have set ourselves the task of determining whether there is more to gravitational 
physics than is revealed by general relativity. Extra dimension theories in general, and 
braneworld models in particular, are an indispensable avenue with which to approach un- 
derstand gravity, post-Einstein. Extra dimensions provide an approach to modifying gravity 
with out abandoning Einstein's general relativity altogether as a basis for understanding the 
fundamental gravitational interaction. Furthermore, the braneworld paradigm (Fig. 1) al- 
lows model builders a tool by which to avoid very real constraints on the number of extra 
dimensions coming from standard model observations. By explicitly pinning matter and 
standard model forces onto a (3-1-1) -dimensional brane Universe while allowing gravity to 
explore the larger, higher-dimensional space, all nongravitational physics follows the stan- 
dard phenomenology. Ultimately the game in braneworld theories is to find a means by 
which to hide the extra dimensions from gravity as well. Gravity is altered in those regimes 
where the extra dimensions manifest themselves. If we wish to explain today's cosmic accel- 
eration as a manifestation of extra dimensions, it makes sense to devise a braneworld theory 
where the extra dimensions are revealed at only the largest of observable distance scales. 
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braneworld 




FIG. 1. DGP gravity employs the braneworld scenario. Matter and all standard model forces 
and particles are pinned to a strictly four-dimensional braneworld. Gravity, however, is free to 
explore the full five-dimensional bulk. 

A. The Formal Arena 

The braneworld theory [4] of Dvali, Gabadadze, and Porrati (DGP) represents a leading 
model for understanding cosmic acceleration as a manifestation of new gravity. The bulk is 
this model is an empty five- dimensional Minkowski space; all energy-momentum is isolated 
on the four-dimensional brane Universe. The theory is described by the action [4]: 

^(5) = -^M' J d'x^ R + J d^x^^) Cm + Sgh . (2.1) 

M is the fundamental five- dimensional Planck scale. The first term in 5(5) is the Einstein- 
Hilbert action in five dimensions for a five-dimensional metric qab (bulk metric) with Ricci 
scalar R and determinant g. The metric g^^} is the induced (four-dimensional) metric on 
the brane, and is its determinant.^ The contribution Sgh to the action is a pure 
divergence necessary to ensure proper boundary conditions in the Euler-Lagrange equations. 
An intrinsic curvature term is added to the brane action [4]: 

-T^M| / d'x4^) . (2.2) 

lOTT J 

Here, Mp is the observed four-dimensional Planck scale. ^ 



^Throughout this paper, we use A, B, . . . = {0, 1, 2, 3, 5} as bulk indices, i^, . . . = {0, 1, 2, 3} as 
brane spacetime indices, and i, j, . . . = {1, 2, 3} as brane spatial indices. 

^Where would such a term come from? The intrinsic curvature term may be viewed as coming from 
effective-action terms induced by quantum matter fluctuations that live exclusively on the brane 
Universe (see [4-6] for details). There is an ongoing discussion as to whether this theory is unstable 
to further quantum gravity corrections on the brane that reveal themselves at phenomenologically 
important scales [7-14]. While several topics covered here are indeed relevant to that discussion 
(particularly Sec. IV), rather than becoming embroiled in this technical issue, we studiously avoid 
quantum gravity issues here and treat gravity as given by Eqs. (2.1) and (2.2) for our discussion, 
and that for cosmological applications, classical gravity physics is sufficient. 
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The gravitational field equations resulting from the action Eqs. (2.1) and (2.2) are 

M^Gab + Ml5 (x^ - z{x^)) G^AB = Stt 5 (x^ - zix^")) Tab{x^) , (2.3) 

where Gab is the five-dimensional Einstein tensor, G^^^ is the Einstein tensor of the induced 
metric on the brane gj^J, and where is the extra spatial coordinate and z{x^) represents 
the location of the brane as a function of the four-dimensional coordinates of our brane 
Universe, {x^}. Note that the energy-momentum tensor only resides on the brane surface, 
as we have constructed. 

While the braneworld paradigm has often been referred to as "string-inspired," we are 
not necessarily wedded to that premise. One can imagine a more conventional scenario where 
physics is still driven by field theory, where the brane is some sort of solitonic domain wall 
and conventional particles and standard model forces are states bound to the domain wall 
using usual quantum arguments. This approach does require that DGP gravity still exhibits 
the same properties as described in this review when the brane has a nonzero thickness 
[15-18,9]. While the situation seems to depend on the specifics of the brane's substrucutre, 
there exist specific scenarios in which it is possible to enjoy the features of DGP gravity 
with a thick, soliton-like brane. 

Unlike other braneworld theories, DGP gravity has a fanciful sort of appeal; it uncannily 
resembles the Flatland-like world one habitually envisions when extra dimensions are in- 
voked. The bulk is large and relatively flat enjoying the properties usually associated with a 
Minkowski spacetime. Bulk observers may look down from on high upon the brane Universe, 
which may be perceived as being an imbedded surface in this larger bulk. It is important to 
note that the brane position remains fully dynamical and is determined by the field equa- 
tions, Eqs. (2.3). While a coordinate system may devised in which the brane appears flat, 
the brane's distortion and motion are, in that situation, registered through metric variables 
that represent the brane's extrinsic curvature. This is a technique used repeatedly in order 
to case the mathematical description of the brane dynamics. Nevertheless, we will often 
refer to a brane in this review as being warped or deformed or the like. This terminology 
is just shorthand for the brane exhibiting a nonzero extrinsic curvature while imagining a 
different coordinate system in which the brane is nontrivially imbedded. 



B. Preliminary Features 

In order to get a qualitative picture of how gravity works for DGP branewolrds, let us take 
small metric fluctuations around flat, empty space and look at gravitational perturbations, 
Hab, where 

Qab = Vab + hAB , (2.4) 
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where tjab is the five-dimensional Minkowski metric. Choosing the harmonic gauge in the 
bulk 

d^hAB = \dBhi , (2.5) 

where the //5-components of this gauge condition leads to /i^5 = so that the surviving 
components are /i^,^ and h^^. The latter component is solved by the following: 

□ (5)/i5 ^ □(5)^/. ^ ^2.6) 

where D^^^ is the five-dimensional d'Alembertian. The //i/-component of the field equations 
Eqs. (2.3) become, after a little manipulation [6], 

M=^n(^) V + M|5(x5)(n(4) V - d^^d^hl) = Stt (t^, - Irj^^T^) S{x') , (2.7) 

where D*^^-* is the four-dimensional (brane) d'Alembertian, and where we take the brane to 
be located at = 0. Fourier transforming just the four- dimensional spacetime to corre- 
sponding momentum coordinates p^, and applying boundary conditions that force gravita- 
tional fluctuations to vanish one approaches spatial infinity, then gravitational fiuctuations 
on the brane take the form [6] 



(2.8) 



+ 2M3p 

We may recover the behavior of the gravitational potentials from this expression. 
There exists a new physical scale, the crossover scale 

Ml , , 

that governs the transition between four-dimensional behavior and five- dimensional behav- 
ior. Ignoring the tensor structure of Eq. (2.8) until future sections, the gravitational potential 
of a source of mass m is 

Fgrav ~ -^^^ , (2.10) 



when r <C tq. When r ro 



V^grav ~ — , (2- 11) 



where the gravitational strengths are given by Gbuik = M"^ and Gbranc = ^p- I-e., the 
potential exhibits four-dimensional behavior at short distances and five-dimensional behav- 
ior (i.e., as if the brane were not there at all) at large distances. For the crossover scale to 
be large, we need a substantial mismatch between Mp, the conventional four-dimensional 
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FIG. 2. At distances much smaller than the crossover scale rg, gravity appears four-dimensional. 
As a graviton propagates, it's amplitude leaks into the bulk. On scales comparable to ro that 
amplitude is attenuated significantly, thus, revealing the extra dimension. 

Planck scale (corresponding to the usual Newton's constant, Gbrane = G) and the funda- 
mental, or bulk, Planck scale M. The fundamental Planck scale M has to be quite small^ 
in order for the energy of gravity fluctuations to be substantially smaller in the bulk versus 
on the brane, the energy of the latter being controlled by Mp. Note that when M is small, 
the corresponding Newton's constant in the bulk, Gbuik? is large. Paradoxically, for a given 
source mass m, gravity is much stronger in the bulk. 

There is a simple intuition for understanding how DGP gravity works and why the 
gravitational potential has its distinctive form. When Mp ^ M, there is a large mismatch 
between the energy of a gravitational fluctuation on the brane versus that in the bulk. I.e., 
imagine a gravitational field of a given amplitude (unitless) and size (measured in distance). 
The corresponding energies are roughly 

^brane ~ M|, ^ d^x{dhf ^ Ml ^ size (2.12) 

Sxdx^idKf ~ X (size) 2 ~ Ebrane X . (2.13) 

What happens then is that while gravitational fluctuations and field are free to explore the 
entire five-dimensional space unfettered, they are much less substantial energetically in the 
bulk. Imagine an analogous situation. Consider the brane as a metal sheet immersed in air. 
The bulk modulus of the metal is much larger than that of air. Now imagine that sound 



^To have tq ~ -^0^^) today's Hubble radius, one needs M ~ 100 MeV. I.e., bulk quantum gravity 
effects will come into play at this low an energy. How does a low-energy quantum gravity not 
have intolerable effects on standard model processes on the brane? Though one does not have a 
complete description of that quantum gravity, one may argue that there is a decoupling mechanism 
that protects brane physics from bulk quantum effects [6]. 
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FIG. 3. A point source of mass m is located on the brane. The gravitational constant on the 
brane is Gbrane = Mp'^, whereas in the bulk Gbuik = (left diagram). Gravitational equipo- 

tential surfaces, however, are not particularly pathological. Near the matter source, those surfaces 
are lens-shaped. At distances farther from the matter source, the equipotential surfaces become 
increasingly spherical, asymptoting to equipotentials of a free point source in five-dimensions. I.e., 
in that limit the brane has no influence (right diagram). 

waves represent gravity. If one strikes the metal plate, the sound wave can propagate freely 
along the metal sheet as well as into the air. However, the energy of the waves in the air is 
so much lower than that in the sheet, the wave in the sheet attenuates very slowly, and the 
wave propagates in the metal sheet virtually as if there were no bulk at all. Only after the 
wave has propagated a long distance is there a substantial amount of attenuation and an 
observer on the sheet can tell an "extra" dimension exists, i.e., that the sound energy must 
have been lost into some unknown region (Fig. 2). Thus at short distances on the sheet, 
sound physics appears lower dimension, but at larger distances corresponding to the distance 
at which substantial attenuation has occurred, sound physics appears higher dimensional. 
In complete accord with this analogy, what results in DGP gravity is a model where gravity 
deviates from conventional Einstein gravity at distances larger than tq. 

While a brane observer is shielded from the presence of the extra dimensions at distance 
scales shorter than the crossover scale, tq. But from the nature of the above analogy, it should 
be clear that the bulk is not particularly shielded from the presence of brane gravitational 
fields. In the bulk, the solution to Eq. (2.7) for a point mass has equipotential surfaces 
as depicted in Fig. 3. From the bulk perspective, the brane looks like a conductor which 
imperfectly repels gravitational potential lines of sources existing away from the brane, and 
one that imperfectly screens the gravitational potential of sources located on the brane [19]. 

III. COSMOLOGY AND GLOBAL STRUCTURE 

Just as gravity appears four-dimensional at short distances and five- dimensional at large 
distances, we expect cosmology to be altered in a corresponding way. Taking the qualitative 
features developed in the last section, we can now construct cosmologies from the field 
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equations Eqs. (2.3). We will find that the cosmology of DGP gravity provides an intriguing 
avenue by which one may explain the contemporary cosmic acceleration as the manifestation 
of an extra dimension, revealed at distances the size of today's Hubble radius. 



A. The Modified Friedmann Equation 

The first work relevant to DGP cosmolgy appeared even before the article by Dvali, 
Gabadadze and Porrati, though these studies were in the context of older braneworld theories 
[20-22] . We follow here the approach of Deffayet [23] who first noted how one recovers a self- 
accelerating solution from the DGP field equations Eqs. (2.3). The general time-dependent 
line element with the isometrics under consideration is of the form 



(3.1) 



where the coordinates we choose are r, the cosmological time; A', the spatial comoving 
coordinates of our observable Universe; and the extra dimension into the bulk. The 
three-dimensional spatial metric is the Kronecker delta, 5ij^ because we focus our atten- 
tion on spatially-fiat cosmologies. While the analysis was originally done for more general 
homogeneous cosmologies, we restrict ourselves here to this observationally viable scenario. 

Recall that all energy-momentum resides on the brane, so that the bulk is empty. The 
field equations Eqs. (2.3) reduce to 
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(3.5) 



in the bulk. Prime denotes differentiation with respect to z and dot denotes differentiation 
with respect to r. We take the brane to be located at 2; = 0. This prescription does not 
restrict brane surface dynamics as it is tantamount to a residual coordinate gauge choice. 
Using a technique first developed in Refs. [24-27], the bulk equations, remarkably, may 
be solved exactly given that the bulk energy-momentum content is so simple (indeed, it is 
empty here). Taking the bulk to be mirror {Z2) symmetric across the brane, we find the 
metric components in Eq. (3.1) are given by [23] 
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N ^iT z - 
a 

A = a =p \z\a (3.6) 
5 = 1, 

where there remains a parameter a(r) that is yet to be determined. Note that a(r) = 
A{t,z — 0) represents the usual scale factor of the four-dimensional cosmology in our brane 
Universe. Note that there are two distinct possible cosmologies associated with the choice of 
sign. They have very different global structures and correspondingly different phenomenolo- 
gies as we will illuminate further. 

First, take the total energy-momentum tensor which includes matter and the cosmolog- 
ical constant on the brane to be 

I brane = 5{z) (Hag {-p, p, p, p, 0) . (3.7) 

In order to determine the scale factor a(r), we need to employ the proper boundary condition 
at the brane. This can be done by taking Eqs. (2.3) and integrating those equations just 
across the brane surface. Then, the boundary conditions at the brane requires 



TV 
A' 
~A 



P (3.8) 



3 

= ^(3p + 2p). (3.9) 

z=0 



Comparing this condition to the bulk solutions, Eqs. (3.6), these junction conditions require 
a constraint on the evolution of a(T). Such an evolution is tantamount to a new set of 
Priedmann equations [23]: 

/f=±^ = ?^p(r). (3.10) 

To 3 

and 

p + 3(p + p)// = , (3.11) 

where we have used the usual Hubble parameter H — a/a. The second of these equations is 
just the usual expression of energy-momentum conservation. The first equation, however, is 
indeed a new Friedmann equation that is a modification of the convential four-dimensional 
Friedmann equation of the standard cosmological model. 

Let us examine Eq. (3.10) more closely. The new contribution from the DGP braneworld 
scenario is the introduction of the term ±H/ro on the left-hand side of the Friedmann 
equation. The choice of sign represent two distinct cosmological phases. Just as gravity 
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FIG. 4. The solid curve depicts Eq. (3.10) while the dotted line represents the conventional 
four-dimensional Priedmann equation. Two cosmological phases clearly emerge for any given spa- 
tially-homogeneous energy-momentum distribution. 

is conventional four-dimensional gravity at short scales and appears five-dimensional at 
large distance scales, so too the Hubble scale, H{t), evolves by the conventional Friedmann 
equation at high Hubble scales but is altered substantially as H(t) approaches r^^. 

Figure 4 depicts the new Friedmann equation. Deffayet [23] first noted that there are 
two distinct cosmological phases. First, there exists the phase in Eq. (3.10) employing the 
upper sign which had already been established in Refs. [20-22], which transitions between 

~ p to ~ p^. We refer to this phase as the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) phase. The other cosmological phase corresponds to the lower sign in Eq. (3.10). 
Here, cosmology at early times again behaves according to a conventional four- dimensional 
Friedmann equation, but at late times asymptotes to a brane self-inflationary phase (the 
asymptotic state was first noted by Shtanov [21]). 

In this latter self-accelerating phase, DGP gravity provides an alternative explanation 
for today's cosmic acceleration. If one were to choose the cosmological phase associated 
with the lower sign in Eq. (3.10), and set the crossover distance scale to be on the order of 
Hq^, where Hq is today's Hubble scale, DGP could account for today's cosmic acceleration 
in terms of the existence of extra dimensions and a modifications of the laws of gravity. The 
Hubble scale, H{t), evolves by the conventional Friedmann equation when the universe is 
young and H{t) is large. As the universe expands and energy density is swept away in this 
expansion, H{t) decreases. When the Hubble scale approaches a critical threshold, H{t) 
stops decreasing and saturates at a constant value, even when the Universe is devoid of 
energy-momentum content. The Universe asymptotes to a deSitter accelerating expansion. 
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B. The Brane Worldsheet in a Minkowski Bulk 



It is instructive to understand the meaning of the two distinct cosmological phases in the 
new Friedmann equation Eq. (3.10). In order to do that, one must acquire a firmer grasp 
on the global structure of the bulk geometry Eqs. (3.6) and its physical meaning [23,28]. 
Starting with Eqs. (3.6) and using the technique developed by Deruelle and Dolezel [29] in a 
more general context, an explicit change of coordinates may be obtained to go to a canonical 
five-dimensional Minkowskian metric 

ds^ = dT^ - {dX^f - {dX^f - {dX^f - {dY')\ (3.12) 

The bulk metric in this cosmological environment is strictly fiat, and whatever nontrivial 
spacetime dynamics we experience on the brane Universe is derived from the particular 
imbedding of a nontrivial brane surface evolving in this trivial bulk geometry. Rather like 
the popular physics picture of our standard cosmology, we may think of our Universe as a 
balloon or some deformable surface expanding and evolving in a larger-dimensional space. 
Here in DGP gravity, that picture is hterally true. 

The coordinate transformation from Eqs. (3.6) to the exphcitly fiat metric Eq. (3.12) is 
given by 

1 r ^ a"^ d fd 

2 J d^ dr \a 

^firi^(^), (3.13) 




2 J d^ dr \a 



where = 6ijX^X^. 

For clarity, we can focus on the early universe of cosmologies of DGP braneworlds to get a 
picture of the global structure, and in particular of the four-dimensional big bang evolution at 
early times. Here, we restrict ourselves to radiation domination (i.e., p — |p) such that, using 
Eqs. (3.10) and (3.11), a(T) = r^/^ when H ^ Vq^ using appropriately normalized time units. 
A more general equation of state does not alter the qualitative picture. Equation (3.10) shows 
that the early cosmological evolution on the brane is independent of the sign choice, though 
the bulk structure will be very different for the two different cosmological phases through 
the persistence of the sign in Eqs. (3.6). The global configuration of the brane worldsheet 
is determined by setting ^ = in the coordinate transformation Eq. (3.13). We get 

T = 

Y^ = ry'{'—-l-T]-^r'/'' (3.14) 
X' - 
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brane worldsheet 




FIG. 5. A schematic representation of the brane worldsheet from an inertial bulk reference 
frame. The bulk time coordinate, T, is the vertical direction, while the other directions represent 
all four spatial bulk coordinates, X* and Y^. The big bang is located along the locus = T, 
while the dotted surface is the future lightcone of the event located at the origin denoted by the 
solid dot. The curves on the brane worldsheet are examples of equal cosmological time, r, curves 
and each is in a plane of constant Y^ + T. Figure from Ref. [28] . 

The locus of points defined by these equations, for all (r. A*), satisfies the relationship 

Y+ = ^j:{Xr + lY^, (3.15) 

where we have defined Y± = \(T ± Y^). Note that if one keeps only the first term, the 
surface defined by Eq. (3.15) would simply be the light cone emerging from the origin at 
(T, X*,y^) = 0. However, the second term ensures that the brane worldsheet is timelike 
except along the y+-axis. Moreover, from Eqs. (3.14), we see that 

Y_ = , (3.16) 

implying that K_ acts as an effective cosmological time coordinate on the brane. The F^-axis 
is a singular locus corresponding to r = 0, or the big bang.^ 

This picture is summarized in Figs. 5 and 6. Taking Y^ as its time coordinate, a bulk 
observer perceives the braneworld as a compact, hyperspherical surface expanding relativis- 
tically from an initial big bang singularity. Figure 5 shows a spacetime diagram representing 



^The big bang singularity when r < oo is just the origin y_ = y+ = X* = and is strictly 
pointlike. The rest of the big bang singularity (i.e., when 1+ > 0) corresponds to the pathological 
case when r = oo. 
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FIG. 6. Taking time-slicings of the spacetime diagram shown in Fig. 5 (and now only suppress- 
ing one spatial variable, rather than two) we see that the big bang starts as a pointlike singularity 
and the brane universe expands as a relativistic shockwave from that origin. Resulting from the 
peculiarities of the coordinate transformation, the big-bang persists eternally as a lightlike singu- 
larity (see Fig. 5) so that for any given time slice, one point on the brane surface is singular and 
is moving at exactly the speed of light relative to a bulk observer. 

this picture where the three dimensional hypersurface of the brane Universe is depicted as a 
circle with each time slice. Note that a bulk observer views the braneworld as spatially com- 
pact, even while a cosmo logical brane observer does not. Simultaneously, a bulk observer 
sees a spatially varying energy density on the brane, whereas a brane observer perceives 
each time slice as spatially homogeneous. Figure 6 depicts the same picture as Fig. 5, but 
with each successive image in the sequence representing a single time slice (as seen by a 
bulk observer). The big bang starts as a strictly pointlike singularity, and the brane surface 
looks like a relativistic shock originating from that point. The expansion factor evolves by 
Eq. (3.10) implying that at early times near the big bang, its expansion is indistinguishable 
from a four-dimensional FLRW big bang and the expansion of the brane bubble decelerates. 
However, as the size of the bubble becomes comparable to tq, the expansion of the bubble 
starts to deviate significantly from four-dimensional FLRW. 

Though the brane cosmological evolution between the FLRW phase and the self- 
accelerating phase is indistinguishable at early times, the bulk metric Eqs. (3.6) for each 
phase is quite distinct. That distinction has a clear geometric interpretation: The FLRW 
phase (upper sign) corresponds to that part of the bulk interior to the brane worldsheet, 
whereas the self-accelerating phase (lower sign) corresponds to bulk exterior to the brane 
worldsheet (see Fig. 7). The full bulk space is two copies of either the interior to the brane 
worldsheet (the FLRW phase) or the exterior (the self-accerating phase), as imposed by 
22^symmetry. Those two copies are then spliced across the brane, so it is clear that the full 
bulk space cannot really be represented as imbedded in a fiat space. 
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FIG. 7. The brane surface at a given time instant as seen from a inertial bulk observer. While 
from a brane observer's point of view (observer b), a constant-time slice of the universe is infinite 
in spatial extent, from a bulk observer's point of view, the brane surface is always compact and 
spheroidal (imagine taking a time slice in Fig. 5). That spheroidal brane surface expands at near 
the speed of light from a bulk observer's point of view. In the FLRW phase, a bulk observer 
exists only inside the expanding brane surface, watching that surface expand away from him/her 
(observer -B+). In the self-accelerating phase, a bulk observer only exists outside the expanding 
brane surface, watch that surface expand towards him/her (observer B-). 

It is clear that the two cosmological phases really are quite distinct, particularly at early 
times when the universe is small. In the FLRW phase, the bulk is the tiny interior of a 
small brane bubble. From a bulk observer's point of view, space is of finite volume, and 
he/she witnesses that bulk space grow as the brane bubble expands away from him/her. The 
intriguing property of this space is that there are shortcuts through the bulk one can take 
from any point on the brane Universe to any other point. Those shortcuts are faster than 
the speed of light on the brane itself, i.e., the speed of a photon stuck on the brane surface 
[28]. In the self-accelerating phase, the bulk is two copies of the infinite volume exterior, 
spliced across the tiny brane bubble. Here a bulk observer witnesses the brane bubble 
rapidly expanding towards him/her, and eventually when the bubble size is comparable to 
the crossover scale tq, the bubble will begin to accelerate towards the observer approaching 
the speed of light. Because of the nature of the bulk space, one cannot take shortcuts through 
the bulk. The fastest way from point A to B in the brane Universe is disappointingly within 
the brane itself. 



C. Luminosity Distances and Other Observational Constraints 

How do we connect our new understanding of cosmology in DGP gravity to the real 
(3-1-1) -dimensional world we actually see? Let us focus our attention on the expansion 
history governed by Eq. (3.10) and ask how one can understand this with an eye toward 
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comparison with existing data. In the dark energy paradigm, one assumes that general 
relativity is still valid and that today's cosmic acceleration is driven by a new smooth, 
cosmological component of sufficient negative pressure (referred to as dark energy) whose 
energy density is given by pde and so that the expansion history of the universe is driven 
by the usual Priedmann equation 

H^^^{Pm + Pde) , (3.17) 

the dark energy has an equation of state w — Pde/ Pde, so that 

PDEir) = p2,sa-^(^+-) , (3.18) 

if w is constant; whereas Pde has more complex time dependence if w varies with redshift. 
Dark energy composed of just a cosmological constant {w = —1) is fully consistent with 
existing observational data; however, both w > —1 and w < —1 remain observationally 
viable possibilities [3]. 

We wish to get a more agile feel for how the modified Friedmann equation of DGP 
gravity, Eq. (3.10), 

ro 3 

behaves in of itself as an evolution equation. We are concerned with the situation where 
the Universe is only populated with pressureless energy-momentum constituents, pm, while 
still accelerating in its expansion at late time. We must then focus on the self-accelerating 
phase (the lower sign) so that 

H^-- = -^Pm{t) . 3.19 

ro 3 

Then the effective dark energy density of the modified Friedmann equation is then 

'-^Pf. = - . (3.20) 
3 ro 

The expansion history of this model and its corresponding luminosity distance redshift re- 
lationship was first studied in Ref. [30,31]. 

By comparing this expression to Eq. (3.18), one can mimic a w-model, albeit with a 
time- varying w. One sees immediately that the effective dark energy density attenuates 
with time in the self-accelerating phase. Employing the intuition devised from Eq. (3.18), 
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this implies that the effective w associated with this effective dark energy must always be 
greater than negative one.^ 

What are the parameters in this model of cosmic expansion? We assume that the universe 
is spatially flat, as we have done consistently throughout this review. Moreover, we assume 
that Hq is given. Then one may define the parameter D,m in the conventional manner such 
that 

nM = ^%{l + zf , (3.21) 

where 

= ^ . (3.22) 



It is imperative to remember that while Pm is the sole energy-momentum component in this 
Universe, spatial flatness does not imply Qm = 1- This identity crucially depends on the 
conventional four-dimensional Friedmann equation. One may introduce a new effective dark 
energy component, flro, where 

1 

to resort to anologous identity: 



^ro = z-EI ^ (3.23) 



l = nM + ^ro ■ (3-24) 

This tactic, or something similar, is often used in the literature. Indeed, one may even 
introduce an effective time-dependent Wes{z) = P^de/p^de- Using Eq. (3.11) and the time 
derivative of Eq. (3.19), 

^eff(^) = -YTcT ■ (3-25) 

Taking Eq. (3.19), we can write the redshift dependence of H{z) in terms of the parameters 
of the model: 



H{z) _ 1 
Ho 2 



(3.26) 



^It must be noted that if one were to go into the FLRW phase, rather than self-accelerating 
phase, and relax the presumption that the cosmological constant be zero (i.e., abandon the notion of 
completely replacing dark energy) , then there exists the intriguing possibility of gracefully achieving 
WeS < — 1 without violating the null-energy condition, without ghost instabilities and without a 
big rip [32-34]. 
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FIG. 8. dl{z)/di^'^^^{z) and for a variety of models. The reference model 

is the best-fit flat ACDM, with 17^^ = 0.27. The dashed curves are for the constant-tf models with 
w) = (0.2, -0.765), (0.27, -0.72), and (0.35, -0.69) from top to bottom. The solid cm^ves are 
for the DGP models with the same as the constant-w curves from top to bottom. 

While it seems that Eq. (3.26) exhibits two independent parameters, and tqHq, Eq. (3.24) 
implies that 

1 



1 



(3.27) 



yielding only a single free parameter in Eq. (3.26). 

The luminosity-distance takes the standard form in spatially flat cosmologies: 

dz 



d 



DGP I 



Hiz) 



(3.28) 



using Eq. (3.26). We can compare this distance with the luminosity distance for a constant- 
w dark-energy model 



dl{z) = {l + z) 



Ho^dz 



^Qfjil + zf + (1 - QI,){1 + ^)3{i+-) ' 



(3.29) 



and in Fig. 8 we compare these two luminosity distances, normalized using the best-fit 
ACDM model for a variety of Q% values. 

What is clear from Fig. 8 is that for all practical purposes, the expansion history of DGP 
self-accelerating cosmologies are indistinguishable from constant-w dark-energy cosmologies. 



18 



They would in fact be identical except for the fact that Wea{z) has a clear and specific 
redshift dependence given by Eq. (3.25). The original analysis done in Ref. [30,31] suggests 
that SNIA data favors an low compared to other independent measurements. Such a 
tendency is typical of models resembling w > —1 dark energy. Supernova data from that 
period implies that the best fit Omega\j is 

= 0.181°:°^ , (3.30) 

at the one-sigma level resulting from chi-squared minimization. Equation (3.24) implies that 
the corresponding best-fit estimation for the crossover scale is 

ro = 1.211°:°^ H^' . (3.31) 

Subsequent work using supernova data [35,36,32,37,33,38^0] refined and generalized these 
results. The most recent supernova results [3] are able to probe the deceleration/acceleration 
transition epoch [40] . Assuming a flat universe, this data suggests a best fit 

= 0.21 , (3.32) 

corresponding to a best-fit crossover scale 

ro = 1.26 H^^ . (3.33) 

Similar results were obtained when relaxing flatness or while applying a gaussian prior on 
the matter density parameter. 

Another pair of interesting possibilities for probing the expansion history of the universe 
is using the angular size of high-z compact radio-sources [41] and using the estimated age 
of high-z objects [42] . Both these constraints are predicated on the premise that the only 
meaningful effect of DGP gravity is the alteration of the expansion history. However, if the 
objects are at high enough redshift, this may be a plausible scenario (see Sec. V). 

Finally, one can combine supernova data with data from the cosmic microwave back- 
ground (CMB). Again, we are presuming the only effect of DGP gravity is to alter the 
expansion history of the universe. While that is mostly likely a safe assumption at the last 
scattering surface (again see Sec. V), there are 0(l)-redshift effects in the CMB, such as the 
late-time integrate Sachs- Wolfe effect, that may be very sensitive to alterations of gravity at 
scales comparable to today's Hubble radius. We pursue such issues later in this review. For 
now, however, we may summarize the findings on the simpler presumption [31]. Supernova 
data favors slightly lower values of Jl^ compared to CMB data for a flat universe. However, 
a concordance model with fl^ = 0.3 provided a good fit to both sets (pre-WMAP CMB 
data) with ~ 140 for the full data set (135 data points) with a best fit crossover scale 
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IV. RECOVERY OF EINSTEIN GRAVITY 



Until now, we have ignored the crucial question of whether adopting DGP gravity yields 
anomalous phenomenology beyond the alteration of cosmic expansion history. If we then 
imagine that today's cosmic acceleration were a manifestation of DGP self-acceleration, the 
naive expectation would be that all anomalous gravitational effects of this theory would 
be safely hidden at distances substantially smaller than today's Hubble radius, Hq^, the 
distance at which the extra dimension is revealed. We will see in this section that this 
appraisal of the observational situation in DGP gravity is too naive. 

DGP gravity represents an infrared modification of general relativity. Such theories 
often have pathologies that render them phenomenologically not viable. These pathologies 
are directly related to van Dam-Veltman-Zakharov discontinuity found in massive gravity 
[43-45]. DGP does not evade such concerns: although gravity in DGP is four-dimensional 
at distances shorter than Tq, it is not four-dimensional Einstein gravity - it is augmented by 
the presence of an ultra-light gravitational scalar, roughly corresponding to the unfettered 
fluctuations of the braneworld on which we live. This extra scalar gravitational interaction 
persists even in the limit where Tq ^ —>■ 0. This is a phenomenological disaster which is only 
averted in a nontrivial and subtle matter [46-49] . Let us first describe the problem in detail 
and then proceed to understanding its resolution. 

A. The van Dan— Veltman— Zakharov Discontinuity 

General relativity is a theory of gravitation that supports a massless graviton with two 
degrees of freedom, i.e., two polarizations. However, if one were to describe gravity with a 
massive tensor field, general covariance is lost and the graviton would possess five degrees 
of freedom. The gravitational potential (represented by the quanitity hfj_i, — g^i, — rjf^i,) 
generated by a static source T/^i, is then given by (in three-dimensional momentum space, 

'(9^) = (^.- (4.1) 

for a massive graviton of mass m around a Minkowski-flat background. While similar in 
form to the gravitational potential in Einstein gravity 

-(g^) = (t,. - Iv,.T:) (4.2) 

it nevertheless has a distinct tensor structure. In the limit of vanishing mass, these five 
degrees of freedom may be decomposed into a massless tensor (the graviton), a massless 
vector (a graviphoton which decouples from any conserved matter source) and a massless 
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scalar. This massless scalar persists as an extra degree of freedom in all regimes of the theory. 
Thus, a massive gravity theory is distinct from Einstein gravity, even in the limit where the 
graviton mass vanishes as one can see when comparing Eqs. (4.1) and (4.2). This discrepancy 
is a formulation of the van Dam-Veltman-Zakharov ( VDVZ) discontinuity [43-45] . 

The most accessible physical consequence of the VDVZ discontinuity is the gravitational 
field of a star or other compact, spherically symmetric source. The ratio of the strength 
of the static (Newtonian) potential to that of the gravitomagnetic potential is different for 
Einstein gravity compared to massive gravity, even in the massless limit. Indeed the ratio 
is altered by a factor of order unity. Thus, such effects as the predicted hght deflection by 
a star would be affected significantly if the graviton had even an infinitesimal mass. 

This discrepancy appears for the gravitational field of any compact object. An even more 
dramatic example of the VDVZ discontinuity occurs for a cosmic string. A cosmic string 
has no static potential in Einstein gravity; however, the same does not hold for a cosmic 
string in massive tensor gravity. One can see why using the potentials Eqs. (4.1) and (4.2). 
The potential between a cosmic string with T^^ — diag(T', — T, 0, 0) and a test particle with 
f^, = diag(2M2, 0, 0, 0) is 

TM 

^massless j ^massive i (^■^) 

where the last expression is taken in the limit m ^ 0. Thus in a massive gravity theory, we 
expect a cosmic string to attract a static test particle, whereas in general relativity, no such 
attraction occurs. The attraction in the massive case can be attributed to the exchange of 
the remnant light scalar mode that comes from the decomposition of the massive graviton 
modes in the massless limit. 

The gravitational potential in DGP gravity, Eq. (2.8), has the same tensor structure as 
that for a massive graviton and pcrturbatively has the same VDVZ problem in the limit that 
the graviton linewidth (effectively Tq^) vanishes. Again, this tensor structure is the result 
of an effective new scalar that may be associated with a branc fluctuation mode, or more 
properly, the fluctuations of the extrinsic curvature of the branc. Because, in this theory, the 
brane is tensionless, its fluctuations represent a very light mode and one may seriously ask 
the question as to whether standard tests of scalar-tensor theories, such as light-deflection 
by the sun, already rule out DGP gravity by wide margins. 

It is an important and relevant question to ask. We are precisely interested in the limit 
when — > for all intents and purposes. We want Tq to be the size of the visible Universe 
today, while all our reliable measurements of gravity are on much smaller scales. However, 
the answers to questions of observational relevance are not straightforward. Even in massive 
gravity, the presence of the VDVZ discontinuity is more subtle than just described. The 
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potential Eq. (4.1) is only derived perturbatively to lowest order in /i^j, or T^,^. Vainshtein 
proposed that this discontinuity does not persist in the fully-nonlinear classical theory [50] . 
However, doubts remain [51] since no self-consistent, fully-nonlinear theory of massive tensor 
gravity exists (see, for example, Ref. [52]). 

If the corrections to Einstein gravity remain large even in limit ro — > oo, the phenomenol- 
ogy of DGP gravity is not viable. The paradox in DGP gravity seems to be that while it is 
clear that a perturbative, VDVZ-like discontinuity occurs in the potential somewhere (i.e., 
Einstein gravity is not recovered at short distances), no such discontinuity appears in the 
cosmological solutions; at high Hubble scales, the theory on the brane appears safely like 
general relativity [46]. What does this mean? What is clear is that the cosmological solu- 
tions at high Hubble scales are extremely nonlinear, and that perhaps, just as Vainshtein 
suggested for massive gravity, nonlinear effects become important in resolving the DGP 
version of the VDVZ discontinuity. 

B. Case Study: Cosmic Strings 

We may ask the question of how nonlinear, nonperturbative effects change the potential 
Eq. (2.8), per se. Indeed, as a stark and straightforward exercise, we may ask the question 
in DGP gravity, does a cosmic string attract a static test particle or not in the hmit? We 
will see that corrections remain small and that the recovery of Einstein gravity is subtle and 
directly analogous to Vainstein's proposal for massive gravity. DGP cosmic strings provided 
the first understanding of how the recovery of Einstein gravity occurs in noncosmological 
solutions [47]. Cosmic strings offer a conceptually clean environment and a geometrically 
appealing picture for how nonperturbative effects drive the loss and recover of the Einstein 
limit in DGP gravity. Once it is understood how the VDVZ issue is resolved in this simpler 
system, understanding it for the Schwarzschild-like solution becomes a straightforward affair. 

1. The Einstein Solution 

Before we attempt to solve the full five-dimensional problem for the cosmic string in 
DGP gravity, it is useful to review the cosmic string solution in four-dimensional Einstein 
gravity [53,54]. For a cosmic string with tension T, the exact metric may be represented by 
the line element: 

ds'' = de - dx" - (1 - 2Gr)-' dr' - r'^dcf)^ . (4.4) 

This represents a flat space with a deficit angle AttGT. If one chooses, one can imagine 
suppressing the x-coordinate and imagining that this analysis is that for a particle in (2+1)- 
dimensional general relativity. Equation (4.4) indicates that there is no Newtonian potential 
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(i.e., the potential between static sources arising from goo) between a cosmic string and a 
static test particle. However, a test particle (massive or massless) suffers an azimuthal 
deflection of AnGT when scattered around the cosmic string, resulting from the deficit angle 
cut from spacetime. Another way of interpreting this defiection effect may be illuminated 
through a different coordinate choice. The line element Eq. (4.4) can be rewritten as 

Again, there is no Newtonian gravitational potential between a cosmic string and a static 
test particle. There is no longer an explicit deficit angle cut from spacetime; however, 
in this coordinate choice, the defiection of a moving test particle results rather from a 
gravitomagnetic force generated by the cosmic string. 

In the weak field limit, one may rewrite Eq. (4.5) as a perturbation around fiat space, 
i-G., g^u, = ri^j, + h^y, as a series in the small parameter GT such that 

^00 = ^xa; = (4.6) 
hyy^h,,^AGT\nr , (4.7) 

where r = \/y^ + z'-^ is the radial distance from the cosmic string. So, interestingly, one does 
recover the logarithmic potentials that are expected for codimension-2 objects like cosmic 
strings in (3+l)-dimensions or point particles in (2+l)-dimensions. They appear, however, 
only in the gravitomagnetic potentials in Einstein gravity, rather than in the gravitoelectric 
(Newtonian) potential. 

2. DGP Cosmic Strings: The Weak-Brane Limit 

We wish to fine the spacetime around a perfectly straight, infinitely thin cosmic string 
with a tension T, located on the surface of our brane Universe (see Fig. 9). Alternatively, 
we can again think of suppressing the coordinate along the string so that we consider the 
spacetime of a point particle, located on a two dimensional brane existing in a (3+1)- 
dimensional bulk. As in the cosmological solution, we assume a mirror, 22^symmetry, 
across the brane surface at 2; = |. The Einstein equations Eqs. (2.3) may now be solved for 
this system. 

There is certainly a regime where one may take a perturbative limit when GT is small 
and so that given g^B = Vab + h^B, the four-dimensional Fourier transform of the metric 
potential on the brane is given by Eq. (2.8). For a cosmic string, this implies that when 
r 3> ro. 
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braneworld 





FIG. 9. A schematic representation of a spatial slice through a cosmic string located at A. The 
coordinate x along the cosmic string is suppressed. The coordinate r represents the 3-dimensional 
distance from the cosmic string A, while the coordinate z denotes the polar angle from the vertical 
axis. In the no-gravity limit, the braneworld is the horizontal plane, z = ^. The coordinate (p is 
the azimuthal coordinate. Figure from Ref. [47]. 
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Graviton modes localized on the brane evaporate into the bulk over distances comparable 
to tq. The presence of the brane becomes increasingly irrelevant as r/rg —> oo and a cosmic 
string on the brane acts as a codimension-three object in the full bulk. When r <^ tq, 
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(4.10) 
(4.11) 



The metric potentials when r <^ ro represent a conical space with deficit angle |47rGT. 
Thus in the weak field limit, we expect not only an extra light scalar field generating the 
Newtonian potential, but also a discrepancy in the deficit angle with respect to the Einstein 
solution. 

We can ask the domain of validity of the perturbative solution. The perturbative solution 
considered only terms in Eqs. (2.3) linear in hAB-, or correspondingly, linear in GT . When 
GT <^ 1, this should be a perfectly valid approach to self-consistenly solving Eqs. (2.3). 
However, there is an important catch. While GT is indeed a small parameter, DGP gravity 
introduces a large parameter tq into the field equations. Actually, since ro is dimensionful, 
the large parameter is more properly rg/r. Thus, there are distances for which nonlinear 
terms in Eqs. (2.3) resulting from contributions from the extrinsic curvature of the brane 



^(GT)2 
r 



(4.12) 



cannot be ignored, even though they are clearly higher order in GT. Nonlinear terms such 
as these may only be ignored when [47] 
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(4.13) 



Thus, the perturbative solution given by the metric potential Eq. (2.8) is not valid in all 
regions. In particular, the perturbative solution is not valid in the limit where everything is 
held fixed and ro — >^ oo, which is precisely the limit of interest. 



3. The r/ro — > Limit 

For values of r violating Eq. (4.13), nonlinear contributions to the Einstein tensor become 
important and the weak field approximation breaks down, even when the components h^i, -C 
1. What happens when r <^ ro\/4jrGT7 We need to find a new small expansion parameter 
in order to find a new solution that applies for small r. Actually, the full field equations 
Eqs. (2.3) provide a clue [47]. A solution that is five-dimensional Ricci fiat in the bulk, 
sporting a brane surface that is four-dimensional Ricci fiat, is clearly a solution. Figure 10 is 
an example of such a solution (almost). The bulk is pure vanilla five-dimensional Minkowski 
space, clearly Ricci fiat. The brane is a conical deficit space, a space whose intrinsic curvature 
is strictly zero. The field equations Eqs. (2.3) should be solved. 

Why the space depicted in Fig. 10 is not exactly a solution comes from the Z2-symmetry 
of the bulk across the brane. The brane surface has nontrivial extrinsic curvature even 
though it has vanishing intrinsic curvature. Thus a polar deficit angle space has a residual 
bulk curvature that is a delta-function at the brane surface, and Eqs. (2.3) are not exactly 
zero everywhere for that space. Fortunately, the residual curvature is subleading in r/ro, 
and one may perform a new systematic perturbation in this new parameter, r/ro, starting 
with the space depicted in Fig. 10 as the zeroth-order limit. 

The new perturbative solution on the brane is given using the line element 

ds'' = N\r)\sim=(3 {dt'' - dx^) - ^'(r)|sin.=/3 dr^ - (5h^d(t>'' , (4.14) 
where the metric components on the brane are [47] 

A^(r)U.=, = l + ^^^f (4.15) 

^(r)U.=, = l-^^^^-, (4.16) 
zp ro 

and the deficit polar angle in the bulk is 7r(l — a) where sin ^ = while the deficit 
azimuthal angle in the brane itself is 27r(l — /3). The deficit angle on the brane is given by 

/3 = 1-2^7, (4.17) 

which is precisely equivalent to the Einstein result. The perturbative scheme is valid when 
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FIG. 10. A spatial slice through the cosmic string located at A. As in Fig. 9 the coordinate x 
along the cosmic string is suppressed. The solid angle wedge exterior to the cone is removed from 
the space, and the upper and lower branches of the cone are identified. This conical surface is the 
braneworld (2; = ^ or sinz = (3). The bulk space now exhibits a deficit polar angle (cf. Fig. 9). 
Note that this deficit in polar angle translates into a conical deficit in the braneworld space. Figure 
from Ref. [47]. 
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FIG. 11. The Newtonian potential for a cosmic string has the following regimes: outside rg, the 
cosmic string appears as a codimension-3 object, i.e., a Schwarzschild source, and so its potential 
goes as r~^; inside ro, the string appears as a codimension-2 object, i.e. a true string source. 
Outside ro(T/Mp)^/^, however, the theory appears Brans-Dicke and one generates a logarithmic 
scalar potential associated with codimension-2 objects. Inside the radius ro(T/Mp)^/^ from the 
string source, Einstein gravity is recovered and there is no Newtonian potential. 



(4.18) 



which is complementary to the regime of validity for the weak-brane perturbation. Moreover, 
Eq. (4.18) is the regime of interest when taking tq — oo while holding everything else fixed, 
i.e., the one of interest in addressing the VDVZ discontinutiy. What we see is that, just like 
for the cosmological solutions, DGP cosmic strings do no suffer a VDVZ-like discontinuity. 
Einstein gravity is recovered in the ro oo limit, precisely because of effects nonlinear, and 
indeed nonperturbative, in GT. 

4- The Picture 



Figure 11 depicts how in different parametric regimes, we find different qualitative be- 
haviors for the brane metric around a cosmic string in DGP gravity. Though we have not 
set out the details here, the different perturbative solutions discussed are part of the same 
solution in the bulk and on the brane [47], i.e., the trunk and the tail of the elephant, as it 
were. For an observer at a distance r ^ tq from the cosmic string, where rQ^ characterizes 
the graviton's effective linewidth, the cosmic string appears as a codimension-three object 
in the full bulk. The metric is Schwarzschild-like in this regime. When r -C ro, brane ef- 
fects become important, and the cosmic string appears as a codimension-two object on the 
brane. If the source is weak (i.e., GT is small), the Einstein solution with a deficit angle 
of AttGT holds on the brane only when r <^ r^yj AttGT . In the region on the brane when 
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r ^ roy/AnGT (but still where r <^ ro), the weak field approximation prevails, the cosmic 
string exhibits a nonvanishing Newtonian potential and space suffers a deficit angle different 
from AnGT. 

The solution presented here supports the Einstein solution near the cosmic string in 
the limit that ro — > oo and recovery of Einstein gravity proceeded precisely as Vainshtein 
suggested it would in the case of massive gravity: nonperturbative effects play a crucial role 
in suppressing the coupling the extra scalar mode. Far from the source, the gravitational 
field is weak, and the geometry of the brane (i.e., its extrinsic curvature with respect to 
the bulk) is not substantially altered by the presence of the cosmic string. The solution is 
a perturbation in GT around the trivial space depicted in Fig. 9. Propagation of the light 
scalar mode is permitted and the solution does not correspond to that from general relativity. 
However near the source, the gravitational fields induce a nonperturbative extrinsic curvature 
in the brane, in a manner reminiscent of the popular science picture used to explain how 
matter sources warp geometry. Here, the picture is literally true. The solution here is a 
perturbation in t/tq around the space depicted in Fig. 10. The brane's extrinsic curvature 
suppresses the coupling of the scalar mode to matter and only the tensor mode remains, 
thus Einstein gravity is recovered. 

C. The Schwarzschild-like Solution 

So while four-dimensional Einstein gravity is recovered in a region near a cosmic string 
source, it is recovered in a region much smaller than the region where one naively expected 
the extra dimension to be hidden, i.e., the larger radius Tq. Einstein gravity is only recov- 
ered within a region much smaller than tq, a region whose size is dictated by the source 
strength. Do the insights elucidated using cosmic string translate for a star-hke object? If 
so, that would have fantastic observational consequences. We would have a strong handle 
for observing this theory in a region that is accessible in principle, i.e., at distances much 
smaller than today's Hubble radius. Indeed, Gruzinov first showed how that recovery of 
Einstein gravity in the Schwarzschild-like solution is exactly analogous to what was found 
for the cosmic string, and moreover, is also in exactly the spirit of Vainshtein's resolution of 
the VDVZ discontinuity for massive gravity [48] . 

1. The Field Equations 

We are interested in finding the metric for a static, compact, spherical source in a 
Minkowski background. Under this circumstance, one can choose a coordinate system in 
which the metric is static (i.e., has a timelike Killing vector) while still respecting the spher- 
ical symmetry of the matter source. Let the line element be 
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ds^ = Ar2(r, z)de - A'(r, z)dr'' - B^{r, z)[d0^ + sin' 0d(f)^] - dz" 



(4.19) 



This is the most general static metric with spherical symmetry on the brane. The bulk 
Einstein tensor for this metric is: 
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(4.20) 



The prime denotes partial differentiation with respect to r, whereas the subscript z represents 
partial differentiation with respect to z. 

We wish to solve the five-dimensional field equations, Eq. (2.3). This implies that all 
components of the Einstein tensor, Eqs. (4.20), vanish in the bulk but satisfy the following 
modified boundary relationships on the brane. Fixing the residual gauge B\z=q — r and 
imposing Z2-symmetry across the brane 
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when z = 0. These brane boundary relations come from Gu, G„- and Ggg, respectively. We 
have chosen a gauge in which the brane, while still dynamical, appears flat. All the important 
extrinsic curvature effects discussed in the last section will appear in the ^-derivatives of the 
metric components evaluated at the brane surface, rather than through any explicit shape 
of the brane. 
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We are interested in a static matter distribution p(r), and we may define an effective 
radially-dependent Schwarzschild radius 
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(4.22) 



where we will also use the true Schwarzschild radius, Vg = Rg{r ^ oo). We are interested 
only in weak matter sources, Pg{r). Moreover, we are most interested in those parts of 
spacetime where deviations of the metric from Minkowski are small. Then, it is convenient 
to define the functions {n(r, z), a{r, z), b{r, z)} such that 



N{r, = 1 + n(r, z) 

A{r, z) = 1 + a{r, z) 
B{r, z) = r [1 + b{r, z)] 



(4.23) 
(4.24) 
(4.25) 



Since we are primarily concerned with the metric on the brane, we can make a gauge choice 
such that 6(r, z — 0) — identically so that on the brane, the hue element 



ds' = [1 + n{r) \,=of - [1 + a{r)\,=oY dr' - r'dVt , 



2 J„2 



(4.26) 



takes the standard form with two potentials, n(r)|^=o ^-nd a{r)\z=o., the Newtonian potential 
and a gravitomagnetic potential. Here we use dVt as shorthand for the usual differential solid 
angle. 

We will be interested in small deviations from flat Minkowski space, or more properly, we 
are only concerned when n(r, 2;), a(r, z) and 6(r, z) <^1. Wc can then rewrite our field equa- 
tions, Eqs. (4.20), and brane boundary conditions, Eqs. (4.21), in terms of these quantities 
and keep only leading orders. The brane boundary conditions become 
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(4.27) 



Covariant conservation of the source on the brane allows one to ascertain the source pressure, 
p(r), given the source density p(r): 

Pg = -n'pg . (4.28) 
The pressure terms were dropped from Eqs. (4.27) because there are subleading here. 
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2. The Weak-Brane Limit 



Just as for the cosmic string, there is again a regime where one may take a properly take 
the perturbative hmit when Vg is smalL Again, given qab = VAB + hAB, the four- dimensional 
Fourier transform of the metric potential on the brane is given by Eq. (2.8) 



For a Schwarzschild solution, this implies that when r ^ ro 
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hxx — hyy — hzz — 'o~2 ' (4.32) 

It is convenient to write the latter in terms of our new potentials for the line element, 
Eq. (4.26), 

^W|.=o = -^J (4.33) 

«(r)|.=o = +^g:. (4.34) 

This is actually the set of potentials one expects from Brans-Dicke scalar-tensor gravity with 
Dicke parameter uj — Einstein gravity would correspond to potentials whose values are 
—rg/2r and -\-rG/2r, respectively. As discussed earher, there is an extra hght scalar mode 
coupled to the matter source. That mode may be interpreted as the fluctuations of the free 
brane surface. 

Again, just as in the cosmic string case, we see that significant deviations from Einstein 
gravity yield nonzero contributions to the right-hand side of Eqs. (4.21). Because these 
are mulitplied by tq, this implies that the extrinsic curvatures (as represented by the z- 
derivatives of the metric components at z = 0) can be quite large. Thus, while we neglected 
nonlinear contributions to the field equations, Eqs. (4.20), bilinear terms in those equations 
of the form for example, are only negligible when [48] 

r:^r,= [rgr^^^''' , (4.35) 

even when Vg is small and n, a <^ 1. When r <^ r*, we need to identify a new perturbation 
scheme. 
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weak-brane phase 




braneworld 



Einstein pliase 



FIG. 12. Given a mass source located on the brane, inside the radius (the green hemisphere), 
the brane is dimpled significantly generating a nonperturbative extrinsic curvature. Brane fluc- 
tuations are suppressed in this region (Einstein phase). Outside this radius, the brane is free to 
fluctuate. The bulk in this picture is above the brane. The mirror copy of the bulk space below 
the brane is suppressed. 



The key to identifying a solution when r <^ r^, is to recognize that one only needs to keep 
certain nonlinear terms in Eqs. (4.20). So long as n,a,b <^ 1, or equivalently r ^ Vg, the 



Consider a point mass source such that Rgir) = Tg = constant. Then, the following set of 
potentials on the brane are [48] 



The full bulk solution and how one arrives at that solution will be spelled out in Sec. V 
when we consider the more general case of the Schwarzschild-like solution in the background 
of a general cosmology, a subset of which is this Minkowski background solution. 

That the inclusion of terms only nonlinear in and hz was sufficient to find solutions 
valid when r <^ r^, is indicative that the nonlinear behavior arises from purely spatial 
geometric factors [48]. In particular, inserting the potentials Eqs. (4.36) and (4.37) into 
the expressions Eqs. (4.27) indicates that the extrinsic curvatures of the brane, i.e., az\z=o 
and hz\z=0i play a crucial role in the nonlinear nature of this solution, indeed a solution 
inherently nonperturbative in the source strength Tg. This again is directly analogous to the 
cosmic string but rather than exhibiting a conical distortion, the brane is now cuspy. The 
picture of what happens physically to the brane is depicted in Fig. 12. When a mass source 
is introduced in the brane, its gravitational effect includes a nonperturbative dimpling of the 
brane surface (in direct analogy with the popular physics picture of how general relativity 
works). The brane is dimpled significantly in a region within a radius r^, of the matter source. 



3. The r/r^ Limit 



only nonlinear terms that need to be included are those terms bilinear in ^ and [48]. 




(4.36) 



(4.37) 
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FIG. 13. The Newtonian potential V(r) = gQQ — 1 has the following regimes: outside rg, 
the potential exhibits five-dimensional behavior (i.e., l/r^); inside tq, the potential is indeed 
four-dimensional (i.e., 1/r) but with a coefficient that depends on r. Outside r^, we have 
Brans-Dicke potential while inside r^, we have a true four-dimensional Einstein potential. 

The extrinsic curvature suppresses the light brane bending mode associated with the extra 
scalar field inside this region, whereas outside this region, the brane bending mode is free 
to propagate. Thus four- dimensional Einstein gravity is recovered close the mass source, 
but at distances less than r*, not distances less than tq. Outside r^,, the theory appears like 
a four-dimensional scalar-tensor theory, in particular, four-dimensional linearized Brans- 
Dicke, with parameter u = 0. A marked departure from Einstein gravity persists down to 
distances much shorter than tq. Figure 13 depicts the hierarchy of scales in this system. 



So, we expect a marked departure of the metric of a spherical, compact mass at distances 
comparable to and greater. The potentials Eqs. (4.36) and (4.37) provide the form of the 
corrections to Einstein gravity as r approaches while the in the weak-brane phase, i.e., 
when r ^ r^, but when r is still much smaller than tq, the potentials are given by Eqs. (4.34) 
and (4.34). From our treatment of the cosmic expansion history of DGP gravity, if we are to 
comprehend the contemporary cosmic acceleration as a manifestation of extra-dimensional 
effects, wish to set tq ~ Hq^. The distance clearly plays an important role in DGP 
phenomenology. Table I gives a few examples of what would be if given source masses 
were isolated in an empty Universe when tq ~ -f^o'^. 

However, there is a complication when we wish to understand the new gravitational forces 
implied by Eqs. (4.34) and (4.34) the context of cosmology [55]. The complication arises 
when we consider cosmologies whose Hubble radii, H~^, are comparable or even smaller 



V. MODIFIED GRAVITATIONAL FORCES 
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TABLE I. Example Values for r* 



Earth 
Sun 

Milky Way (10 




1.2 pc 
150 pc 
1.2 Mpc 



than Tq. Take H ^ = constant = Tq for example. In such an example, one may actually 
regard the Hubble flow as being as making an effective contribution to the metric potentials 



Figure 14 depicts a representative situation. The corrections computed in the Minkowski case 
become important just at the length scales where the cosmology dominates the gravitational 
potential. I.e., the regime r ^ r* is just that region where cosmology is more important 
than the localized source. Inside this radius, an observer is bound in the gravity well of the 
central matter source. Outside this radius, an observer is swept away into the cosmological 
flow. Thus, one cannot reliably apply the results from a Minkowski background under 
the circumstance of a nontrivial cosmology, particularly, when we are interested in DGP 
gravity because of its anomalous cosmological evolution. We need to redo the computation 
to include the background cosmology, and indeed we will find a cosmology dependent new 
gravitational force [55,56,34]. 

This computation is more nuanced than for a static matter source in a static Minkowski 
background. We are still interested in finding the metric for compact, spherically symmetric 
overdensities, but now there is time-evolution to contend with. However, if we restrict 
our attention to distance scales such that rH <^ 1 and to nonrelativistic matter sources,^ 
then to leading-order in r^if^ and zH, the solutions to the field equations Eqs. (2.3) are 
also solutions to the static equations, i.e. the metric is quasistatic, where the only time 
dependence comes from the slow evolution of the extrinsic curvature of the brane. To be 
explicit, we are looking at the nonrelativistic limit, where the gravitational potentials of a 
matter source depend only on the instantaneous location of the matter's constituents, and 
not on the motion of those constituents. Incidentally, we left the details of the solution to 
the Minkowski problem (which were first treated in Refs. [48,49]) to this section. All the 
arguments to be employed may be used, a fortiori, as a special case of the more general 
cosmological background. 



^There are a number of simplifications that result from these approximations. See Ref. [56] for 
an enumeration of these as well as the caveats concerning straying to far from the approximations. 



cosmo — 



(5.1) 
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FIG. 14. The corrections to the Newtonian potential become as large as the potential itself 
just as the contributions from cosmology become dominant. It is for this reason that we expect 
the background cosmology to has a significant effect on the modified potential in the weak-brane 
regime. 



35 



A. Background Cosmology 



One can choose a coordinate system in which the cosmological metric respects the spher- 
ical symmetry of the matter source. We are concerned with processes at distances, r, such 
that rH ^ 1. Under that circumstance it is useful to change coordinates to a frame that 
surrenders explicit brane spatial homogeneity but preserves isotropy 

r(r,A^) =a(r)A (5.2) 
t{T,y)=T+^H{T)a'{T) , (5.3) 

for all z and where — SijXX^ . The hne element becomes 



1^2{H + H/H)\z\ - {H^ + Hy] dt^ - [1 ^ 2H\z\] [(1 + H\^)dr^ + r^dU] - dz^ , 

(5.4) 



where here dot repreresents differentiation with respect to the new time coordinate, t. More- 
over, H = H{t) in this coordinate system. All terms of 0{r^H^) or 0{z'^H'^, zHrH) and 
higher have been neglected. The key is that because we are interested primarily in phenom- 
ena whose size is much smaller than the cosmic horizon, the effect of cosmology is almost 
exclusively to control the extrinsic curvature, of the brane. This can he interpreted as a 
modulation of the brane 's stiffness or the strength of the scalar gravitational mode. 

In the coordinate system described by Eq. (5.4), the bulk is like a Rindlcr space. This has 
a fairly natural interpretation if one imagines the bulk picture [23,28]. One imagines riding 
a local patch of the brane, which appears as hyper spherical surface expanding into (or away 
from) a five-dimensional Minkowski bulk. This surface cither accelerates or decelerates in its 
motion with respect to the bulk, creating a Rindler-typc potential. Note that we arc keeping 
the upper-and-lower-sign convention to represent the two cosmological phases. While we are 
nominally focussed on self-acceleration, we will see that contributions from the sign have 
important effects on the modified gravitational potentials. 



B. Metric Potentials 



We have chosen a coordinate system, Eq. (4.19), in which a compact spherical matter 
source may have a quasistatic metric, yet still exist within a background cosmology that is 
nontrivial (i.e., deSitter expansion). Let us treat the matter distribution to be that required 
for the background cosmology, Eq. (5.4), and add to that a compact spherically symmetric 
matter source, located on the brane around the origin {r — 0, z — 0) 

r^lbrane = S{z) diag {6p{r) + pb, -5p{r) - pb, -Sp{r) - pb, -5p{r) - pb, 0) , (5.5) 
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where ps and ps are the density and pressure of the background cosmology, and where 
5p{r) is the overdensity of interest and 5p{r) is chosen to ensure the matter distribution and 
metric are quasistatic. We may define an effective Schwarzschild radius 



(5.6) 



We solve the perturbed Einstein equations in quasistatic approximation by generalizing the 
method used in [55], obtaining the metric of a spherical mass overdensity Sp{t,r) in the 
background of the cosmology described by Eqs. (3.10) and (3.11). Because we are interested 
only in weak matter sources, Pg{r), and since we are interested in solutions well away from 
the cosmic horizon, we can still exploit doing an expansion using Eqs. (4.24) and keeping 
only leading orders. Now, we just need to take care to include leading orders in terms of H 
as well as other parameters of interest. 

We are particularly concerned with the evaluation of the metric on the brane. However, 
we need to take care that when such an evaluation is performed, proper boundary conditions 
in the bulk are satisfied, i.e., that there no singularities in the bulk, which is tantamount 
to having spurious mass sources there. In order to determine the metric on the brane, we 
implement the following approximation [55,56,34]: 



(5.7) 



Note that this is just the value Uz would take if there were only a background cosmology, 
but we are making an assumption that the presence of the mass source only contributes 
negligibly to this quantity at the brane surface. With this one specification, a complete set 
of equations, represented by the brane boundary conditions Eqs. (4.21) and = 0, exists 
on the brane so that the metric functions may be solved on that surface without reference 
to the bulk. At the end of this section, we check that the assumption Eq. (5.7) indeed turns 
out the be the correct one that ensures that no pathologies arise in the bulk. 
The brane boundary conditions Eqs. (4.21) now take the form 
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where we have substituted Eqs. (3.10) and (3.11) for the background cosmological density 
and pressure and where we have neglected second-order contributions (including those from 
the pressure necessary to keep the compact matter source quasistatic). 
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FIG. 15. The function A(r) represents a normalized correction to Newton's constant, G, i.e., 
Geff = G[l + A(r)]. In the self-accelerating cosmological phase, for small r, A(r) asymptotes to 
— (r^/2rgrQ)^/^, i.e., a correction independent of cosmology. For large r (but also when r <C H~^), 
A(r) asymptotes the constant value 1/3/3. This value is — ^ in the saturated limit vqH = 1, and 
goes like 0{l/rQH) as roH — > cxd. The boundary between the two regimes is thus r^, = (rgrg//?^)^/^. 
For the FLRW phase, the graph is just changed by a sign flip, with the exception that the most 
extreme curve occurs not when rgi/ = 1, but rather when H = 0. 
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There are now five equations on the brane with five unknowns. The solution on the brane 
is given by the following. For a cosmological background with arbitrary evolution H{t), we 
find that [55,56,34] 



rn'it, r)|brane = |^ [1 + A(r)] - {H' + Hy 
zr 

a(t,r)|brane = §[l-A(r)] + ii/V^ 



(5.9) 
(5.10) 



Note that the cosmological background contribution is included in these metric components. 
The function A(r) is defined as 



A(r) 
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\ 9/52r3 



(5.11) 



and 



p^l± 2roH 1 + 



3IP 



(5.12) 



Just as for the modified Friedmann equation, Eq. (3.10), there is a sign degeneracy in 
Eq. (5.12). The lower sign corresponds to the self-accelerating cosmologies. These expres- 
sions are valid on the brane when r ^ rg, H~^. In both Eqs. (5.9) and (5.10), the first term 
represent the usual Schwarzschild contribution with a correction governed by A(r) result- 
ing from brane dynamics (as depicted in Fig. 15), whereas the second term represents the 
leading cosmological contribution. Let us try to understand the character of the corrections. 



C. Gravitational Regimes 



One may consolidate all our results and show from Eqs. (5.9)-(5.11) that there exists a 
scale [55,56,34], 



1/3 



(5.13) 



with 



/3 = 1 ± 2roH 1 + 



JL 
3IP 



Inside a radius r* the metric is dominated by Einstein gravity but has corrections which 
depend on the global cosmological phase. Outside this radius (but at distances much smaller 
than both the crossover scale, ro, and the cosmological horizon, H^^) the metric is weak- 
brane and resembles a scalar-tensor gravity in the background of a deSitter expansion. This 
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braneworld 




FIG. 16. The four-dimensional universe where we hve is denoted by the large spherical brane. 
A local mass source located, for example, near its north pole dynamically dimples the brane, 
inducing a nonperturbative extrinsic curvature. That extrinsic curvature suppress coupling of the 
mass source to the extra scalar mode and, within the region dictated by the radius r=K given by 
Eq. (5.13), Einstein gravity is recovered. Outside r^, the gravitational field is still modulated by 
the effects of the extrinsic curvature of the brane generated by the background cosmology. 

scale is modulated both by the nonperturbative extrinsic curvature effects of the source itself 
as well as the extrinsic curvature of the brane generated by its cosmology. The qualitative 
picture described in Fig. 12 is generalized to the picture shown in Fig. 16. 

Keeping Fig. 15 in mind, there are important asymptotic limits of physical relevance for 
the metric on the brane, Eqs. (5.9) and (5.10). First, when r -C r^,, the metric is close to the 
Schwarzschild solution of four- dimensional general relativity. Corrections to that solution 
are small: 



The background cosmological expansion becomes largely unimportant and the corrections 
are dominated by effects already seen in the Minkowski background. Indeed, there is no 
explicit dependence on the parameter governing cosmological expansion, H. However, the 
sign of the correction to the Schwarzschild solution is dependent on the global properties of 
the cosmological phase. Thus, we may ascertain information about bulk, five- dimensional 
cosmological behavior from testing details of the metric where naively one would not expect 
cosmological contributions to be important. 




(5.14) 



(5.15) 
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Cosmological effects become important wfien r ^ r*. Tfie metric is dominated by tfie 
cosmological flow, but tliere is still an attractive potential associated with the central mass 
source. Within the cosmological horizon, r ^ H~^, this residual potential is 



2r 
2r 



1 

~ 3^ 



(5.16) 
(5.17) 



This is the direct analog of the weak-brane phase one finds for compact sources in Minkowski 
space. The residual effect of the matter source is a linearized scalar-tensor gravity with 
Brans-Dicke parameter 

3. 



(5.18) 



Notice that as roi? — > oo, we recover the Einstein solution, corroborating results found 
for linearized cosmological perturbations [62]. At redshifts of order unity, r^H ~ 0(1) and 
corrections to four-dimensional Einstein gravity are substantial and i7(2;)-dependent. These 
results found in this analysis were also followed up in Ref . [63] . 

D. Bulk Solutions 



We should elaborate how the solution Eqs. (5.9) and (5.10) were ascertained as well as 
explicitly write the solutions of the potentials {n, a, 6} in the bulk [55,56,34]. The key to our 
solution is the assumption Eq. (5.7). It allows one to solve the brane equations independently 
of the rest of the bulk and guarantees the asymptotic bulk boundary conditions. In order 
to see why Eq. (5.7) is a reasonable approximation, we need to explore the full solution to 
the bulk Einstein equations, 

G'As(r,^) = 0, (5.19) 

satisfying the brane boundary conditions Eqs. (5.8), as well as specifying that the metric 
approach the cosmological background Eq. (5.4) for large values of r and i.e., far away 
from the compact matter source. 

First, it is convenient to consider not only the components of the Einstein tensor 
Eqs. (4.20), but also the following components of the bulk Ricci tensor (which also van- 
ishes in the bulk) : 
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(5.20) 
(5.21) 



A B 

We wish to take Gzr = 0, = 0, and i?^ = and derive expressions for A{r, z) and B{r, z) 
in terms of N{r,z). Only two of these three equations are independent, but it is useful to 
use all three to ascertain the desired expressions. 
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1. Weak-Field Analysis 



Since we are only interested in metric when r, 2; -C ro, for a weak matter source, 
we may rewrite the necessary field equations using the expressions Eqs. (4.24). Since the 
functions, {n(r, z), a{r. z), h{r, z)} are small, we need only keep nonlinear terms that include 
^-derivatives. The brane boundary conditions, Eqs. (4.21), suggest that and terms 
may be sufficiently large to warrant inclusion of their nonlinear contributions. It is these z- 
derivative nonlinear terms that are crucial to the recover of Einstein gravity near the matter 
source. If one neglected these bilinear terms as well, one would revert to the linearized, 
weak-brane solution (cf. Ref. [48]). 

Integrating Eq. (5.21) twice with respect to the 2;-coordinate, we get 

n + a + 2h = zgi{r) + g2{r) , (5.22) 

where gi{r) and g2{T) are to be specified by the brane boundary conditions, Eqs. (5.8), and 
the residual gauge freedom 5h{r)\z=Q — 0, respectively. Integrating the Gx^-component of 
the bulk Einstein tensor Eqs. (4.20) with respect to the 2;-coordinate yields 

r (n + 26)' - 2 (a - 6) = g^{r) . (5.23) 

The functions gi{r), 5*2 ('^), and (73 (r) are not all independent, and one can ascertain their re- 
lationship with one another by substituting Eqs. (5.22) and (5.23) into the bulk equation. 
If one can approximate = ^{H + H/H) for all 2;, then one can see that Gzr = 0, = 0, 
and Rl = are all consistently satisfied by Eqs. (5.22) and (5.23), where the functions gi{r), 
g2{r), and g3{r) are determined at the brane using Eqs. (5.9) and (5.10) and the residual 
gauge freedom b{r)\z=o ~ 0: 

^1^ = ^(4//+ I) -^(7?, A)' (5.24) 

g2{r) = -^H + ^(1 - A) + c^r ^(1 + A) (5.25) 

g^ir) = ^(1 - 3A) - {2H^ + i?) r^ (5.26) 

where we have used the function A(r), defined in Eq. (5.11). Using Eqs. (5.22)-(5.26), we 
now have expressions for a(r, z) and 6(r, z) completely in terms of n(r, z) for all (r, z). 

Now we wish to find n(r, z) and to confirm that = =F(if + H/H) is a good approxi- 
mation everywhere of interest. Equation (5.20) becomes 



„ 2n' ,H^ + H 

n H h rizz = ±- 



- H 



9i{r) ± — — 



(5.27) 
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where again we have neglected contributions if we are only concerned with r, z <^ ro, H ^. 
Using the expression Eq. (5.24), we find 



n" + — + n,, = -3 {h' + h)t 2^ . (5.28) 



Then, if we let 



n 



It (^H + z - ^(h^ + h) T ro ^ dr ^Rg{r)A{r) + 6n{r, z) , (5.29) 



where 5n{r, z) satisfies the equation 



5n" + ^ + 5n,, = , (5.30) 
r 



we can solve Eq. (5.30) by requiring that 5n vanish as r, ^ ^ oo and applying the condition 

(5.31) 



r on L-o = 



l+(^l±2ro^^^) A(r) 



on the brane as an alternative to the appropriate brane boundary condition for Sn{r, z) 
coming from a linear combination of Eqs. (5.8). We can write the solution explicitly: 

5n(r,z)^ dk cik)e-^' sinkr , (5.32) 
Jo 

where 

2 

c{k) = — J dr r sin kr Sn\ (r) . (5.33) 
We can then compute SuzI^^q, arriving at the bound 

Sn.\z-o^- Tdr^^ , (5.34) 
r Jo H 

for all r < Tq, H-\ Then, 

n,|,^o = IF + + 5n,|,^o . (5.35) 

When the first term in Eq. (5.35) is much larger than the second, Eq. (5.7) is a good 
approximation. When the two terms in Eq. (5.35) are comparable or when the second term 
is much larger than the first, neither term is important in the determination of Eqs. (5.9) 
and (5.10). Thus, Eq. (5.7) is still a safe approximation. 

One can confirm that all the components of the five- dimensional Einstein tensor, 
Eqs. (4.20), vanish in the bulk using field variables satisfying the relationships Eqs. (5.22), 
(5.23), and (5.29). The field variables a{r,z) and b{r,z) both have terms that grow with 
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stemming from the presence of the matter source. However, one can see that with the 
following redefinition of coordinates: 

R = r- ztq-^ (5.36) 

Z^z+Tdr^, (5.37) 
Jo 

that to leading order as 2; — > H~^, the desired Z-dependence is recovered for a{R, Z) and 
b{R, Z) (i.e., ^HZ), and the Newtonian potential takes the form 

n{R, Z)^^[h + ^Z -]^{H^ + H)R^ + ■■■ . (5.38) 

Thus, we recover the desired asymptotic form for the metric of a static, compact matter 
source in the background of a cosmological expansion. 



2. A Note on Bulk Boundary Conditions 



A number of studies have been performed for DGP gravity that have either arrived 
at or used modified force laws different than those given by Eqs. (5.9)-(5.12) [57-61]. How 
does one understand or reconcile such a discrepancy? Remember that one may ascertain the 
metric on the brane without reference to the bulk because there are five unknown quantities, 
{iV(r), A{r), Nz{r), Az{r), i?^(r)}, and there are four independent equations (the three brane 
boundary conditions, Eqs. (5.8), and the G^-component of the bulk Einstein equations) on 
the brane in terms of only these quantities. One need only choose an additional relationship 
between the five quantities in order to form a closed, solvable system. We chose Eq. (5.7) 
and showed here that it was equivalent to choosing the bulk boundary condition that as 
z became large, one recovers the background cosmology. The choice made in these other 
analyses is tantamount to a different choice of bulk boundary conditions. One must be 
very careful in the analysis of DGP gravitational fields that one is properly treating the 
asymptotic bulk space, as it has a significant effect on the form of the metric on the brane. 



VI. ANOMALOUS ORBIT PRECESSION 



We have established that DGP gravity is capable of generating a contemporary cosmic 
acceleration with dark energy. However, it is of utmost interest to understand how one 
may differentiate such a radical new theory from a more conventional dark energy model 
concocted to mimic a cosmic expansion history identical to that of DGP gravity. The 
results of the previous sections involving the nontrivial recovery of Einstein gravity and 
large deviations of this theory from Einstein gravity in observably accessible regimes is the 
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key to this observational differentiation. Again, there are two regimes where one can expect 
to test this theory and thus, there are two clear regimes in which to observationally challenge 
the theory. First deep within the gravity well, where r <S r* and where the corrections to 
general relativity are small, but the uncertainties are also correspondingly well-controlled. 
The second regime is out in the cosmological flow, where r ^ r^^ (but still r ^ ro ~ H^^) 
and where corrections to general relativity are large, but our observations are also not as 
precise. We focus on the first possibility in this section and will go to the latter in the next 
section. 



A. Nearly Circular Orbits 

Deep in the gravity well of a matter source, where the effects of cosmology are ostensibly 
irrelevant, the correction to the gravitational potentials may be represented by effective 
correction to Newton's constant 



A(r) = ±, 



(6.1) 



that appears in the expressions Eqs. (5.9) and (5.10). Though there is no explicit dependence 
on the background H{t) evolution, there is a residual dependence on the cosmological phase 
through the overall sign. In DGP gravity, tests of a source's Newtonian force leads to 
discrepancies with general relativity [48,55,64]. 

Imagine a body orbiting a mass source where Rg{r) — Vg — constant. The perihelion pre- 
cession per orbit may be determined in the usual way given a metric of the form Eqs. (4.24) 
and (4.26) 

[■ J AN , , 

^<l>- dr - , (6.2) 

J ^ ^£;2 _ (i + ^) 

where E — N'^dt/ds and J = r'^d(l)/ds are constants of motion resulting from the isometrics of 
the metric, and ds is the differential proper time of the orbiting body. With a nearly circular 
orbit deep within the Einstein regime (i.e., when r <^ r* so that we may use Eqs. (5.14) and 
(5.15)), the above expression yields 

A0 = 2.+ -^^-^^j . (6.3) 

The second term is the famous precession correction from general relativity. The last term 
is the new anomalous precession due to DGP brane effects. This latter correction is a purely 
Newtonian effect. Recall that any deviations of a Newtonian central potential from 1/r 
results in an orbit precession. The DPG correction to the precession rate is now [48,55,64] 
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d 3 

^^0DGP = = T5 /xas/year . (6.4) 

Note that this result is independent of the source mass, implying that this precession rate is a 
universal quantity dependent only on the graviton's effective linewidth (vq^) and the overall 
cosmological phase. Moreover, the final result depends on the sign of the cosmological phase 
[55] . Thus one can tell by the sign of the precession whether self- acceleration is the correct 
driving force for today's cosmic acceleration. It is extraordinary that a local measurement, 
e.g. in the inner solar system, can have something definitive to say about gravity on the 
largest observable scales. 



B. Solar System Tests 



Nordtvedt [65] quotes precision for perihelion precession at 430 //as/ycar for Mercury and 
10 //as/year for Mars. Improvements in lunar ranging measurements [66,64] suggest that the 
Moon will be sensitive to the DGP correction Eq. (6.4) in the future lunar ranging studies. 
Also, BepiColombo (an ESA sateUite) and MESSENGER (NASA) are being sent to Mercury 
at the end of the decade, will also be sensitive to this correction [67]. Counterintuitively, 
futTirc outer solar system probes that possess precision ranging instrumentation, such as 
Cassini [68] , may also provide ideal tests of the anomalous precession. Unlike post-Newtonian 
deviations arising from Einstein corrections, this anomaly does not attenuate with distance 
from the sun; indeed, it amplifies. This is not surprising since we know that corrections 
to Einstein gravity grow as gravity weakens. More work is needed to ascertain whether 
important inner or outer solar system systematics allow this anomalous precession effect to 
manifest itself effectively. If so, we may enjoy precision tests of general relativity even in 
unanticipated regimes [69-73] . The solar system seems to provide a most promising means 
to constrain this anomalous precession from DGP gravity.^ 

it is also interesting to contrast solar system numbers with those for binary pulsars. 
The rate of periastron advance for the object PSR 1913+16 is known to a precision of 
4 X 10'' /xas/year [75]. This precision is not as good as that for the inner solar system. A 
tightly bound system such as a binary pulsar is better at finding general relativity corrections 
to Newtonian gravity because it is a stronger gravity system than the solar system. It is 
for precisely the same reason that binary pulsars prove to be a worse environment to test 



^One should also note that DGP gravity cannot explain the famous Pioneer anomaly [74]. First, 
the functional form of the anomalous acceleration in DGP is not correct: it decays as a r~^/^ 
power-law with distance from the sun. Moreover, the magnitude of the effect is far too small in 
the outer solar system. 
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^^^^^^ cosmic horizon of A 




FIG. 17. The universe is populated by a variety of matter sources. For sources massive enough 
to deform the brane substantially, Einstein gravity is recovered within the gravity well of each 
region (the green circles whose sizes are governed by ~ (mass)^/^). Outside these gravity wells, 
extra scalar forces play a role. 

DGP effects, as these latter effects become less prominent as the gravitational strength of 
the sources become stronger. 

As a final note, one must be extremely careful about the application of Eqs. (5.9) and 
(5.10). Remember they were derived for a spherical source in isolation. We found that the 
resulting metric, while in a weak- field regime, i.e., r^/r -C 1, was nevertheless still very 
nonlinear. Thus, superposibility of central sources is no longer possible. This introduces 
a major problem for the practical application of Eqs. (5.9) and (5.10) to real systems. If 
we take the intuition developed in the last two sections, however, we can develop a sensible 
picture of when these equations may be applicable. We found that being in the Einstein 
regime corresponds to being deep within the gravity well of a given matter source. Within 
that gravity well, the extrinsic curvature of the brane generated by the source suppressed 
coupling of the extra scalar mode. Outside a source's gravity well, one is caught up in 
the cosmological flow, the brane is free to fluctuate, and the gravitational attraction is 
augmented by an extra scalar field. 

If one takes a "rubber sheet" picture of how this DGP effect works, we can imagine 
several mass sources placed on the brane sheet, each deforming the brane in a region whose 
size corresponds with the mass of the source (see Fig. 17). These deformed regions will in 
general overlap with each other in some nontrivial fashion. However, the DGP corrections 
for the orbit of a particular test object should be dominated by the gravity well in which the 
test object is orbiting. For example, within the solar system, we are clearly in the gravity 
well of the sun, even though we also exist inside a larger galaxy, which in turn is inside 
a cluster, etc. Nevertheless, the spacetime curvature around us is dominated by the sun, 
and for DGP that implies the extrinsic curvature of the brane in our solar system is also 
dominated by the sun. This picture also implies that if we are not rooted in the gravity 



47 



well of a single body, then the quantitative form of the correction given by Eqs. (5.14) and 
(5.15) is simply invalid. I.e., three body systems need to be completely reanalyzed. This 
may have relevant consequences for the moon which has substantial gravitational influences 
from both the earth and the sun. 

VII. LARGE-SCALE STRUCTURE 

Future solar system tests have the possibility of probing the residual deviation from four- 
dimensional Einstein gravity at distances well below r*. Nevertheless, it would be ideal to 
test gravitational physics where dramatic differences from Einstein gravity are anticipated. 
Again, this is the crucial element to the program of differentiating a modified-gravity scenario 
such as DGP from a dark-energy explanation of today's cosmic acceleration that has an 
identical expansion history. A detailed study of large scale structure in the Universe can 
provide a test of gravitational physics at large distance scales where we expect anomalous 
effects from DGP corrections to be large. In the last section, we have seen that the modified 
force law is sensitive to the background cosmological expansion, since this expansion is 
intimately tied to the extrinsic curvature of the brane [23,28], and this curvature controls 
the effective Newtonian potential [55,56,34]. This gives us some measure of sensitivity to 
how cosmology affects the changes in the growth of structure through the modulation of 
a cosmology-dependent Newtonian potential. We may then proceed and compare those 
results to the standard cosmology, as well as to a cosmology that exactly mimics the DGP 
expansion history using dark energy. The description of the work presented in this section 
first appeared in Ref. [56]. 

A. Linear Growth 

The modified gravitational potentials given by Eqs. (5.9) and (5.10) indicate that sub- 
stantial deviations from Einstein gravity occur when r ^ r*. More generally, from our 
discussion of the VDVZ discontinuity, we expect large deviations from general relativity any 
time an analysis perturbative in the source strength is valid. This is true when we want 
to understand the linear growth of structure in our universe. So long as we consider early 
cosmological times, when perturbations around the homogeneous cosmology are small, or at 
later times but also on scales much larger than the clustering scale, a linear analysis should 
be safe. 

The linear regime is precisely what we have analyzed with the potential given by Eq. (2.8), 
or, more generally, Eqs. (5.16) and (5.17) 
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in a background cosmology (remember, we are still only considering nonrelativistic sources 
at scales smaller than the horizon size). Because we are in the linear regime, we may think 
of these results as Green's function solutions to the general physical problem. Thus, these 
potentials are applicable beyond just spherical sources and we deduce that 



V^(5n(r, t) = 47rG 



1 



Sp{x, t) 



(7.1) 
(7.2) 



for general matter distributions, so long as they arc nonrelativistic. What results is a 
linearized scalar-tensor theory with cosmology- dependent Brans-Dicke parameter Eq. (5.18) 



Note this is only true for weak, linear perturbations around the cosmological background, 
and that once overdensities becomes massive enough for self-binding, these results cease to 
be valid and one needs to go to the nonlinear treatment. 

Nonrelativistic matter perturbations with S{t) — Sp/p{t) evolve via 



6 + 2Hd = A-nGp 



(7.3) 



where p{t) is the background cosmological energy density, implying that self-attraction of 
overdensities is governed by an evolving G^s'- 



GeS — G 



(7.4) 



Here the modification manifests itself as a time-dependent effective Newton's constant, G^s- 
Again, as we are focused on the self-accelerating phase, then from Eq. (5.12) 



/3 = 1 - 2roH 1 



3IP 



As time evolves the effective gravitational constant decreases. For example, if = 0.3, 
Geff/G = 0.72, 0.86, 0.92 at ^ = 0, 1, 2. 

One may best observe this anomalous repulsion (compared to general relativity) through 
the growth of large-scale structure in the early universe. That growth is governed not 
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FIG. 18. The top panel shows the ratio of the growth factors (dashed hnes) in DGP gravity 
[Eq. (7.3)] and a model of dark energy (DE) with an equation of state such that it gives rise to the 
same expansion history (i.e. given by Eq. (3.10), but where the force law is still given by general 
relativity). The upper line corresponds to = 0.3, the lower one to = 0.2. The sohd lines 
show the analogous result for velocity perturbations factors / = din D-^-/ din a. The bottom panel 
shows the growth factors as a function of redshift for models with different expansion histories, 
corresponding to (from top to bottom) ACDM {Q^ = 0.3), and DGP gravity with = 0.3,0.2 
respectively. Figure from Ref. [56]. 




only by the expansion of the universe, but also by the gravitational self-attraction of small 
overdensities. Figure 18 depicts how growth of large-scale structure is altered by DGP 
gravity. The results make two important points: 1) growth is suppressed compared to the 
standard cosmological model since the expansion history is equivalent to a w{z) dark-energy 
model with an effective equation of state given by Eq. (3.25) 

1 



WeffUj 



and 2) growth is suppressed even compared to a dark energy model that has identical 
expansion history from the alteration of the self attraction from the modified Newton's 
constant, Eq. (7.4). The latter point reiterates the crucial feature that one can differentiate 
between this modified-gravity model and dark energy. 
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B. Nonlinear Growth 



Wc certainly want to understand large scale structure beyond the linear analysis. Unlike 
the standard cosmological scenario where the self-gravitation of overdensities may be treated 
with a linear Newtonian gravitational field, in DGP gravity the gravitational field is highly 
nonlinear, even though the field is weak. This nonlinearity, particularly the inability to con- 
vert the gravitational field into a superposition of point-to-point forces, poses an enormous 
challenge to understanding growth of structure in DGP gravity. 

Nevertheless, we already have the tools to offer at least a primitive preliminary analysis. 
We do understand the full nonlinear gravitational field around spherical, nonrelativistic 
sources, Eqs. (5.9)-(5.12). Consider the evolution of a spherical top-hat perturbation S{t,r) 
of top-hat radius Rt- At subhorizon scales {Hr -C 1), the contribution from the Newtonian 
potential, n(t,r), dominates the geodesic evolution of the overdensity. The equation of 
motion for the perturbation is [56] 



where e = ^r^Rg / 9 13"^ . Note that for large 5, Eq. (7.5) reduces to the standard evolution 
of spherical perturbations in general relativity. 

Figure 19 shows an example of a full solution of Eq. (7.5) and the corresponding solution 
in the cosmological constant case. Whereas such a perturbation collapses in the ACDM 
case dX z = 0.66 when its linearly extrapolated density contrast is 5c = 1.689, for the 
DGP case the collapse happens much later at z = 0.35 when its Sc = 1.656. In terms of 
the linearly extrapolated density contrasts things do not look very different, in fact, when 
the full solutions are expressed as a function of the linearly extrapolated density contrasts, 
5iin = -D+5j/(-D+).j they are very similar to within a few percent. This implies that all the 
higher-order moments of the density field are very close to those for ACDM models, for 
example, the skewness is less than a 1% difference from ACDM. 

This close correspondence of the higher-order moments can be useful by allowing the use 
of non-linear growth to constrain the bias between galaxies and dark matter in the same way 
as it is done in standard case, thus inferring the linear growth factor from the normalization 
of the power spectrum in the linear regime. Although the result in the right panel in Fig. 19 
may seem a coincidence at first sight, Eq. (7.5) says that the nontrivial correction from 
DGP gravity in square brackets is maximum when 6 = (which gives the renormalization of 
Newton's constant). As 6 increases the correction disappears (since DGP becomes Einstein 
at high-densities), so most of the difference between the two evolutions happens in the linear 
regime, which is encoded in the linear growth factor. 




(7.5) 
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FIG. 19. Numerical solution of the spherical collapse. The left panel shows the evolution for a 
spherical perturbation with 5j = 3 x 10~^ at Zi = 1000 for = 0.3 in DGP gravity and in ACDM. 
The right panel shows the ratio of the solutions once they are both expressed as a function of their 
linear density contrasts. Figure from Ref. [56]. 

C. Observational Consequences 



What are the implications of these results for testing DGP gravity using large-scale struc- 
ture? A clear signature of DGP gravity is the suppressed (compared to ACDM) growth of 
perturbations in the linear regime due to the different expansion history and the addition 
of a repulsive contribution to the force law. However, in order to predict the present nor- 
malization of the power spectrum at large scales, we need to know the normalization of the 
power spectrum at early times from the CMB. A fit of the to pre-WMAP CMB data was 
performed in Ref. [31] using the angular diameter distance for DGP gravity, finding a best 
fit (fiat) model with ^ 0.3, with a very similar CMB power spectrum to the standard 
cosmological constant model (with ~ 0.3 and Vl\ = 0.7) and other parameters kept fixed 
at the same value. Here we use this fact, plus the normalization obtained from the best- 
fit cosmological constant power-law model from WMAP [76] which has basically the same 
(relevant for large-scale structure) parameters as in Ref. [31], except for the normalization 
of the primordial fluctuations which has increased compared to pre-WMAP data (see e.g. 
Fig. 11 in Ref. [77]). The normalization for the cosmological constant scale-invariant model 
corresponds to present rms fluctuations in spheres of 8 Mpc/h, dg = 0.9 ±0.1 (see Table 2 
in [76]). 
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FIG. 20. The linear power spectrum normalization, erg, for DGP gravity as a function of O^. 
The vertical lines denote the best fit value and 68% confidence level error bars from fitting to 
type-IA supernovae data from [31], = 0.18lo;gg. The other lines correspond to as as a function 
of obtained by evolving the primordial spectrum as determined by WMAP by the DGP growth 
factor. Figure from Ref. [56]. 

Figure 20 shows the present value of erg as a function of for DGP gravity, where we 
assume that the best-fit normaUzation of the primordial fluctuations stays constant as we 
change Q^, and recompute the transfer function and growth factor as we move away from 
fi^ = 0.3. Since most of the contribution to (Xg comes from scales r < 100/i/Mpc, we can 
calculate the transfer function using Einstein gravity, since these modes entered the Hubble 
radius at redshifts high enough that they evolve in the standard fashion. The value of ag at 

= 0.3 is then given by 0.9 times the ratio of the DGP to ACDM growth factors shown in 
the bottom panel of Fig. 18. The error bars in ag reflect the uncertainty in the normalization 
of primordial fluctuations, and we keep them a constant fraction as we vary away from 
0.3. We see in Fig. 20 that for the lower values of preferred by fltting the acceleration 
of the universe, the additional suppression of growth plus the change in the shape of the 
density power spectrum drive ag to a rather small value. This could in part be ameliorated 
by increasing the Hubble constant, but not to the extent needed to keep erg at reasonable 
values. The vertical lines show the best-flt and la error bars from fltting DGP gravity to 
the supernova data from Ref. [31]. This shows that fltting the acceleration of the universe 
requires approximately cxg < 0.7 to la and as < 0.8 to 2a. 
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In order to compare this prediction of erg to observations one must be careful since most 
determinations of as have built in the assumption of Einstein gravity or ACDM models. We 
use galaxy clustering, which in view of the results in Sect. VII B for higher-order moments, 
should provide a test of galaxy biasing independent of gravity being DGP or Einstein. Recent 
determinations of erg from galaxy clustering in the SDSS survey [78] give — 0.89 ± 0.02 
for L* galaxies at an effective redshift of the survey Zg — 0.1. We can convert this value 
to as for dark matter at 2; = as follows. We evolve to z = using a conservative growth 
factor, that of DGP for = 0.2. In order to convert from L* galaxies to dark matter, we 
use the results of the bispectrum analysis of the 2dF survey [79] where b = 1.04 ±0.11 for 
luminosity L ~ 1.9L*. Wc then scale to L* galaxies using the empirical relative bias relation 
obtained in [80] that 6/6* = 0.85 + 0.15(L/L*), which is in very good agreement with SDSS 
(see Fig. 30 in Ref. [78]). This implies erg = 1.00 ± 0.11. Even if we allow for another 10% 
systematic uncertainty in this procedure, the preferred value of in DGP gravity that 
fits the supernovae data is about 2a away from that required by the growth of structure at 
z = 0. 

Nevertheless, the main difficulty for DGP gravity to simultaneously explain cosmic accel- 
eration and the growth of structure is easy to understand: the expansion history is already 
significantly different from a cosmological constant, corresponding to an effective equation 
of state with w^s = — (1 + ^m)^^- This larger value of w suppresses the growth somewhat 
due to earlier epoch of the onset of acceleration. In addition, the new repulsive contribution 
to the force law suppreses the growth even more, driving ag to a rather low value, in contrast 
with observations. If as error bars shrink the supernovae results continue to be consistent 
with Wed = — 1, this will drive the DGP fit to a yet lower value of and thus a smaller 
value of (Tg. For these reasons wc expect the tension between explaining acceleration and the 
growth of structure to be robust to a more complete treatment of the comparison of DGP 
gravity against observations. Already more ideas for how to approach testing Newton's 
constant on ultralarge scale and the self-consistency of the DGP paradigm for explaining 
cosmic accleration have been taken [81-83,61,84]. 

VIII. GRAVITATIONAL LENSING 
A. Localized Sources 

One clear path to differentiating DGP gravity from conventional Einstein gravity is 
through an anomalous mismatch between the mass of a compact object as computed by 
lensing measurements, versus the mass of an object as computed using some measure of the 
Newtonian potential, such as using the orbit of some satellite object, or other means such 
as the source's x-ray temperature or through the SZ-effect. 
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FIG. 21. Mass discrepancy, 5M, for a static point source whose Schwarzschild radius is rg. The 
sohd curve is for a self- accelerating background with H = Tq^. The dashed curve is for a FLRW 
background with H = r^^. Figure from Ref. [55]. 

The lensing of light by a compact matter source with metric Eq. (5.9) and (5.10) may 
be computed in the usual way. The angle of deflection of a massless test particle is given by 



where E = N'^dt/dX and J = r'^d(j)/dX are constants of motion resulting from the isometries, 
and dX is the differential affine parameter. Removing the effect of the background cosmology 
and just focussing on the deflection generated by passing close to a matter source, the angle 
of deflection is [55] 



where b is the impact parameter. This result is equivalent to the Einstein result, implying 
light deflection is unaltered by DGP corrections, even when those corrections are large. 
This result jibes with the picture that DGP corrections come solely from a light scalar mode 
associated with brane fluctuations. Since scalars do not couple to photons, the trajectory 
of light in a DGP gravitational field should be identical to that in a Einstein gravitational 
field generated by the same mass distribution. 

Light deflection measurements thus probe the "true" mass of a given matter distribution. 
Contrast this with a Newtonian measurement of a matter distribution using the gravitational 




(8.1) 
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potential Eq. (5.9) while incorrectly assuming general relativity holds. The mass discrepancy 
between the lensing mass (the actual mass) and that determined from the Newtonian force 
may be read directly from Eq. (5.9), 

SM = - 1 = ^ - 1 . (8.3) 

This ratio is depicted in Fig. 21 for both cosmological phases, with an arbitrary background 
cosmology chosen. When the mass is measured deep within the Einstein regime, the mass 
discrepancy simplifies to 

Solar system measurements are too coarse to be able to resolve the DGP discrepancy between 
lensing mass of the sun and its Newtonian mass. The discrepancy 6M for the sun at 0(AU) 
scale distances is approximately 10~^^. Limits on this discrepancy for the solar system as 
characterized by the post- Newtonian parameter, 7 — 1, are only constrained to be < 3 x 10~^. 

In the solar system, this is much too small an effect to be a serious means of testing 
this theory, even with the most recent data from Cassini [68] . Indeed most state-of-the-art 
constraints on the post-Newtonian parameter 7 will not do, because most of the tests of 
general relativity focus on the strong-gravity regime near r ~ r^. However, in this theory, 
this is where the anomalous corrections to general relativity are the weakest. 

The best place to test this theory is as close to the weak-brane regime as we can get. 
A possibly more promising regime may be found in galaxy clusters. For 10^^ — > 10^^ Mq 
clusters, the scale {rlRgY^^ has the range 6 — > 14 Mpc. For masses measured at the 
cluster virial radii of roughly 1 — > 3 Mpc, this implies mass discrepancies of roughly 5 
S%. X-ray or Sunyaev-Zeldovich (SZ) measurements are poised to map the Newtonian 
potential of the galaxy clusters, whereas weak lensing measurements can directly measure 
the cluster mass profile. Unfortunately, these measurements are far from achieving the 
desired precisions. If one can extend mass measurements to distances on the order of Tq, 
Fig. 21 suggests discrepancies can be as large as —10% for the FLRW phase or even 50% 
for the self-accelerating phase; however, remember these asymptotic limits get smaller as 
{roH)~^ as a function of the redshift of the lensing mass. 

It is a nontrivial result that light deflection by a compact spherical source is identical to 
that in four-dimensional Einstein gravity (even with potentials Eqs. (5.9)-(5.12) substan- 
tially differing from those of Einstein gravity) through the nonlinear transition between the 
Einstein phase and the weak-brane phase. As such, there remains the possibility that for 
aspherical lenses that this surprising null result does not persist through that transition and 
that DGP may manifest itself through some anomalous lensing feature. However, given the 
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intuition that the entirety of the anomalous DGP effect comes from an extra gravitational 
scalar suggests that even in a more general situation, photons should still probe the "true" 
matter distribution. Were photon geodesies in DGP gravity to act just as in general rela- 
tivity for any distribution of lenses, this would provide a powerful tool for proceeding with 
the analysis of weak lensing in DGP. Weak lensing would then provide a clean window into 
the true evolution of structure as an important method of differentiating it from general 
relativity. Tanaka's analysis of gravitational fields in more general matter distributions in 
DGP gravity provides a first step to answering the nature of photon geodesies [85] . 



The latc-timc integrated Sachs-Wolfe (ISW) effect on the cosmic microwave background 
(CMB) may be viewed as another possible "lensing" signature. The Sachs-Wolfe effect 
is intrinsically connected to how photons move through and arc altered by gravitational 
potentials (in this case, timc-dcpcndcnt ones). This effect is a direct probe of any possible 
alteration of the gravitational potentials over time; moreover, it is a probe of potential 
perturbations on the largest of distance scales, safely in the linear regime. 

Late-time ISW would then seem to be a promising candidate for modified-gravity theories 
of the sort that anticipate cosmic acceleration, as we see that the potentials are altered 
by substantial corrections. Indeed, for a particular set of modified-gravity theories, this 
assertion is indeed the case [86]. However, for DGP gravity, there is an unfortunate catch. 

Recall in Sec. VII, we argued that when considering linear potentials satisfied Eqs. (7.1) 
and (7.2), our results generalized beyond just spherical perturbations. In order to bring our 
linear potentials Eqs. (7.1) and (7.2) into line with how the late-time ISW is normally treated, 
we need to identify the gravitational potentials as perturbations around a homogeneous 
cosmological background with the line element 



Here $(t. A) and ^(i. A) are the relevant gravitational potentials and A is a comoving radial 
coordinate. In effect we want to determine $ and ^ given 5n and 5a. Unhke the case of 
Einstein's gravity, $ ^ One may perform a coordinate transformation to determine 

that relationship. We find that, assigning r = a(T)A, and 



B. The Late-Time ISW 



ds^ = [1 + 2$(t, A)] dr^ - a2(r) [1 + 2^(t, A)] + X^d^ 



(8.5) 



$ = 5n(r, r) 



(8.6) 
(8.7) 




keeping only the important terms when rH -C 1. 
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The quantity of interest for the ISW effect is the time derivative of $ — the above the 
above analysis imphes 

V^($-^) = ^aV, (8.8) 

where V is the gradient in comoving spatial coordiantcs, A*. Just as we found for light 
deflection, this result is identical to the four-dimcnsional Einstein result, the contributions 
from the brane effects exactly cancelling. Again, the intuition of the anomalous DGP effects 
coming from a light gravitational scalar is correct in suggesting the microwave background 
photons probe the "true" matter fluctuations on the largest of scales. 

Thus, the late-time ISW effect for DGP gravity will be identical to that of a dark energy 
cosmology that mimics the DGP cosmic expansion history, Eq. (3.10), at least at scales 
small compared to the horizon. Our approximation does not allow us to address the ISW 
effect at the largest scales (relevant for the CMB at low multipoles), but it is applicable 
to the cross-correlation of the CMB with galaxy surveys [87,88]. At larger scales, one 
expects to encounter difficulties associated with leakage of gravity off the brane (for order- 
unity redshifts) and other bulk effects [28,62,89] that we were successfully able to ignore at 
subhorizon scales. 



C. Leakage and Depletion of Anisotropic Power 

There is an important effect we have completely ignored up until now. At scales com- 
parable to To, when the Hubble expansion rate is comparable to r^^, i.e., 0{1) redshifts, 
gravitational perturbations can substantially leak off the brane. This was the original effect 
discussed from the introduction of DGP braneworlds. At the same time, perturbations that 
exist in the bulk have an equal likelihood of impinging substantial perturbation amplitude 
onto the brane from the outside bulk world. 

This leads to a whole new arena of possibilities and headaches. There is httle that can be 
done to control what exists outside in the free bulk. However, there are possible reasonable 
avenues one can take to simplify the situation. In Sec. IIIB we saw how null worldlines 
through the bulk in the FLRW phase could connect different events on the brane. This 
observation was a consequence of the convexity of the brane worldsheet and the choice of 
bulk. Conversely, if one chooses the bulk corresponding to the self- accelerating phase, one 
may conclude that no null lightray through the bulk connects two different events on the 
brane. 

Gravity in the bulk is just empty space five- dimensional Einstein gravity, and thus, 
perturbations in the bulk must follow null geodesies. Consider again Fig. 5. If we live 
in the self- accelerating cosmological phase, then the bulk exists only exterior to the brane 
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worldsheet in this picture. One can see that, unhke in the interior phase, null geodesies 
can only intersect oTir brane Universe once. I.e., once a perturbation has left our brane, 
it can never get back. Therefore, if we assume that the bulk is completely empty and all 
perturbations originate from the brane Universe, that means that perturbations can only 
leak from the brane, the net result of which is a systematic loss of power with regard to the 
gravitational fluctuations. 

Let us attempt to quantify this depletion. As a crude estimate, we can take the propaga- 
tor from the linear theory given by Eqs. (2.3). and treat only fluctuations in the Newtonian 
potential, 0(a;'^) where g^o = 1 + 20. For modes inside the horizon, we may approximate 
evolution assuming a Minkowski background. For a mode whose spatial wavenumber is k, 
the equations of motion are 



subject to the boundary condition at the brane 

(90 _ _ /a20 



+ = , 



0) 



(8.9) 



(8.10) 



Then for a mode initially localized near the brane, the amplitude on the brane obeys the 
following form 



loe V u , (8.11) 

when T ^ kr^. Imagine the late universe as each mode reenters the horizon. Modes, being 
frozen outside the horizon, are now free to evolve. For a given mode, k, the time spent inside 
the horizon is 

.kro, 

in a late-time, matter-dominated universe and where we have approximated today's cosmic 
time to be tq. Then, the anomalous depletion resulting in DGP gravity is 

1 



1 - 



.12) 







= exp 




DGP 






1 - 



- 1 . 



.13) 



{kro)\ 

This depletion is concentrated at scales where kro ~ 1. It should appear as an enhancement 
of the late-time integrated Sachs- Wolfe effect at the largest of angular scales. 

A more complete analysis of the perturbations of the full metric is needed to get a better 
handle on leakage on sizes comparable to the crossover scale, Tq. Moreover, complicated 
global structure in the bulk renders the situation even more baffling. For example, the 
inclusion of an initial inflationary stage completely alters the bulk boundary conditions. 
Other subtleties of bulk initial and boundary condition also need to be taken into account 
for a proper treatment of leakage in a cosmological setting [28,62,89]. 
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IX. PROSPECTS AND COMPLICATIONS 



We have presented a preliminary assessment of the observational viability of DGP grav- 
ity to simultaneously explain the acceleration of the universe and offer a promising set of 
observables for the near future. The theory is poised on the threshold of a more mature 
analysis that requires a new level of computational sophistication. In order to improve the 
comparison against observations a number of issues need to be resolved. The applicability of 
coarse graining is a pertinent issue that needs to be addressed and tackled. Understanding 
growth of structure into clustered objects is essential in order to apply promising observing 
techniques in the near future. To do a full comparison of the CMB power spectrum against 
data, it remains to properly treat the ISW effect at scales comparable to the horizon. Un- 
derstanding how primordial fluctuations were born in this theory likewise requires a much 
more detailed treatment of technical issues. 

A. Spherical Symmetry and Birkhoff's Law 

An important outstanding question in DGP cosmology is whether a universe driven 
by a uniform distribution of compact sources, such as galaxies or galaxy clusters, actually 
drives the same overall expansion history as a truly uniform distribution of pressureless 
matter. The problem is very close to the same problem as in general relativity (although in 
DGP, the system is more nonlinear) of how one recovers the expansion history of a uniform 
distribution of matter with a statistically homogeneous and isotropic distribution of discrete 
sources. What protects this coarse-graining in general relativity is that the theory possesses 
a Birkhoff's law property (Fig. 22), even in the fully- nonlinear case. 

It is clear from the force law generated from the metric components Eqs. (5.9-5.12) that 
Birkhoff's law does not strictly apply in DGP gravity, even with a pressureless background 
energy-momentum component. The expressions in these equations make a clear distinction 
between the role of matter that is part of an overdensity and matter that is part of the back- 
ground. However, there is a limited sense in which Birkhoff's law does apply. We do see that 
the evolution of overdensities is only dependent on the quantity Rg{r). In this sense, because 
for spherically-symmetric matter configurations, only the integrated mass of the overdensity 
matters (as opposed to more complicated details of that overdensity), Birkhoff's law does 
apply in a limited fashion. So, if all matter perturbations were spherically-symmetric, then 
coarse graining would apply in DGP gravity. 

Of course, matter perturbations are not spherically symmetric, not even when considering 
general relativity. We extrapolate that because we are concerned with perturbations that 
are at least statistical isotropic, that like in general relativity, coarse graining may be applied 
in this more general circumstance. It is widely believed that coarse graining may be applied 
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FIG. 22. It is the BirkhofF's law property of general relativity that allows one take an arbitrary, 
though spherically symmetric matter distribution (left) but describe the evolution of observer B 
relatively to the origin of symmetry A by excising all matter outside a spherical ball of mass 
around the A and condensing the matter within the ball to a single point. DGP gravity, in a 
limited fashion, obeys the same property. 

with impunity in general relativity. Based on the arguments used here, we assume that the 
same holds true for DGP. However, it is worthwhile to confirm this in some more systematic 
and thorough way, as there are subtleties in DGP gravity that do not appear in general 
relativity. 

B. Beyond Isolated Spherical Perturbations 

Understanding the growth of structure in DGP gravity is an essential requirement if one 
is to treat it as a realistic or sophisticated counterpart to the standard cosmological model. 
We need to understand large scale structure beyond a linear analysis and beyond the prim- 
itive spherical collapse picture. Without that understanding, analysis of anything involving 
clustered structures are suspect, particularly those that have formed during the epoch of 
(9(1) redshift where the nonlinear effects are prominent and irreducible. This includes galaxy 
cluster counting, gravitational lens distributions, and several favored methods for determin- 
ing cTg observationally. For analyses of DGP phenomenology such as Refs. [90-96,59,60] to 
carry weight, they need to include such effects. 

What is needed is the moral equivalent of the N-body simulation for DGP gravity. But 
the native nonlinearity coursing throughout this model essentially forbids a simple N-body 
approach in DGP gravity. The "rubber sheet" picture described at the end of Sec. VI and in 
Fig. 17 must be developed in a more precise and formal manner. The truly daunting aspect 
of this approach is that it represents a full five-dimensional boundary value problem for a 
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complex distribution of matter sources, even if one could use the quasistatic approximation. 
One possible simplification is the equivalent of employing the assumption Eq. (5.7) but in a 
more general context of an arbitrary matter distribution. This allows one to reduce the five- 
dimensional boundary value problem to a much simpler (though still terribly intractable) 
four-dimensional boundary value problem. Tanaka has provided a framework where one 
can begin this program of extending our understanding of gravitational fields of interesting 
matter distributions beyond just spherically- symmetric ones [85]. 

If this program can be carried out properly, whole new vistas of approaches can be taken 
for constraining DGP gravity and then this modified-gravity theory may take a place as a 
legitimate alternative to the standard cosmological model. 



There are a number of intriguing exotic possibilities for the future in DGP gravity. 
The model has many rich features that invite us to explore them. One important avenue 
that begs to be developed is the proper inclusion of infiation in DGP gravity. Key issues 
involve the proper treatment of the bulk geometry, how (and at what scale) infiationary 
perturbation in this theory develop and what role do brane boundary conditions play in 
the treatment of this system. An intriguing possibility in DGP gravity may come into 
play as the fundamental five-dimensional Planck scale is so low (corresponding to M ~ 
100 MeV, suggesting transplanckian effects may be important at quite low energies. Initial 
work investigating inflation in DGP gravity has begun [97-99] . Other work that has shown 
just a sample of the richness of DGP gravity and its observational possibilities include the 
possibihty of nonperturbatively dressed or screened starlike solutions [100,101] and shock 
solutions [102,103]. 

Indeed, a possible glimpse of how strange DGP phenomenology may be is exhibited by 
the starlike solution posed in Refs. [100,101]. In these papers, the solution appears as four- 
dimensional Einstein when r <^ r^, just as the one described in Sec. IV. However, far from 
the matter source when r ^ tq, rather than the brane having no influence on the metric, in 
this new solution, there is a strong backrcaction from the brane curvature that screens the 
native mass, M, of the central object and the effective ADM mass appears to be 



where Vg is the Schwarzschild radius corresponding to the mass M. Given the strongly 
nonlinear nature of this system, that two completely disconnected metric solutions exist for 
the same matter source is an unavoidable logical possibihty.^ The consequences of such an 
exotic phenomenology are still to be fully revealed. 



^There is a potential subtlety that may shed light on the nature of this system. A coordinate 



C. Exotic Phenomenology 




) 



1/3 



(9.1) 
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D. Ghosts and Instability 



There has been a suggestion that the scalar, brane-fluctuation mode is a ghost in the 
self- accelerating phase [7,12]. Ghosts are troublesome because, having a negative kinetic 
term, high momentum modes are extremely unstable leading to the rapid devolution of 
the vacuum into increasingly short wavelength fluctuations. This would imply that the 
vacuum is unstable to unacceptable fluctuations, rendering moot all of the work done in 
DGP phenomenology. The status of this situation is still not resolved. Koyama suggests 
that, in fact, around the self-accelerating empty space solution itself, the brane fluctuation 
mode is not a ghost [104], implying that there may be a subtlety involving the self-consistent 
treatment of perturbations around the background. Nevertheless, the situation may be 
even more complicated, particularly in a system more complex than a quasistatic deSitter 
background. 

However, recall that the coupling of the scalar mode to matter is increasingly suppressed 
at sort scales (or correspondingly, high momentum). Even if the extra scalar mode were 
a ghost, it is not clear that the normal mechanism for the instability of the vacuum is as 
rapid here, possibly leaving a phenomenological loophole: with a normal ghost coupled to 
matter, the vacuum would almost instantaneously dissolve into high momentum matter and 
ghosts, here because of the suppressed coupling at high momentum, that devolution may 
proceed much more slowly. A more quantitative analysis is necessary to see if this is a viable 
scenario. 

X. GRAVITY'S FUTURE 

What is gravity? Our understanding of this fundamental, universal interaction has been 
a pillar of modern science for centuries. Physics has primarily focused on gravity as a 
theory incomplete in the ultraviolet regime, at incomprehensibly large energies above the 
Planck scale. However, there still remains a tantalizing regime where we may challenge 
Einstein's theory of gravitation on the largest of scales, a new infrared frontier. This new 
century's golden age of cosmology offers an unprecedented opportunity to understand new 
infrared fundamental physics. And while particle cosmology is the celebration of the intimate 



transformation may take the solution given in Refs. [100,101] into a metric of the form given by 
Eq. (4.19). From this point of view, it is clear that the form of the induced metric on the brane is 
constrained. It would seem then that the imbedding of the brane in the bulk is not fully dynamical, 
and would thus imply that the solution given is not the result of purely gravitational effects. I.e., 
that there must be some extra-gravitational force constraining the brane to hold a particular form. 
A consensus has not yet been achieved on this last conculsion. 



63 



connection between the very great and the very small, phenomena on immense scales in of 
themselves may be an indispensable route to understanding the true nature of our Universe. 

The braneworld theory of Dvali, Gabadadze and Porrati (DGP) has pioneered this new- 
line of investigation, offering an novel explanation for today's cosmic acceleration as resulting 
from the unveiling of an unseen extra dimension at distances comparable today's Hubble 
radius. While the theory offers a specific prediction for the detailed expansion history of the 
universe, which may be tested observationally, it offers a paradigm for nature truly distinct 
from dark energy explanations of cosmic acceleration, even those that perfectly mimic the 
same expansion history. DGP braneworld theory alters the gravitational interaction itself, 
yielding unexpected phenomenological handles beyond just expansion history. Tests from 
the solar system, large scale structure, lensing all offer a window into understanding the 
perplexing nature of the cosmic acceleration and, perhaps, of gravity itself. 

Understanding the complete nature of the gravitational interaction is the final frontier of 
theoretical physics and cosmology, and provides an indispensable and tantalizing opportunity 
to peel back the curtain of the unknown. Regardless of the ultimate explanation, revealing 
the structure of physics at the largest of scales allows us to peer into its most fundamental 
structure. What is gravity? It is not a new question, but it is a good time to ask. 
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